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ABSTRACT 


The  Standard  Ship  Motion  Program,  5MP,  wa3  developed 
at  the  David  Taylor  Naval  Ship  Research  and  Development 
Center  in  1981  as  a  prediction  tool  for  use  in  the  Navy's 
ship  design  process.  SMP  provides  predictions  of  th 
responses  of  a  ship  advancing  at  constant  forward  speed 
with  arbitrary  heading  in  both  regular  waves  and  irregular 
seas.  Since  1981  a  number  of  corrections  and  improvements 
were  made  to  SMP  which  are  detailed  in  this  report. 

Corrections  were  made  in  the  bilge  keel  viscous  damping 
calculation  for  ships  that  have  both  a  bilge  keel  and  a 
skeg  described  on  the  same  station.  An  improvement  in 
yaw-roll  coupling  was  made  by  using  a  different  theory  to 
compute  hull  lift  damping.  New  predictions  were  incorpo¬ 
rated  into  SMP  to  compute  stabilized  ship  responses  for 
active  antiroll  fins  using  fixed  gains.  Predictions  of 
vertical  wave  induced  loads  were  also  incorporated  into  SMP. 

In  addition,  a  new  set  of  tables  of  severe  ship  responses 
is  provided  as  part  of  the  output. 

ADMINISTRATIVE  INFORMATION 

The  David  W.  Taylor  Naval  Ship  Research  and  Development  Center  (DTNSRDC)  was 
authorized  and  funded  over  a  number  of  years  to  develop  and  improve  the  capability 
and  accuracy  of  a  standard  Navy  ship  notion  prediction  computer  program  and  other 
associated  computer  programs.  This  report  documents  computer  program  changes  and 
additions  to  the  report  DTNSRDC /SPD-09 36-01  entitled  "User's  Manual  for  the 
Standard  Ship  Motion  Program,  SMP"  known  as  SMP81,  Funding  was  provided  by  the 
Surface  Ship  Hydromechanics  Program  under  Project  Element  625**3N»  Block 
SF-U3-U21-OOI ,  identified  as  Work  Unit  Numbers  1506-103  and  1506-153.  Funding  was 
also  provided  for  wave  induced  vertical  loads  by  Work  Request  N00Q2483WR1U526 , 
identified  as  Work  Unit  Number  1522-950. 


INTRODUCTION 

The  Standard  Ship  Motion  Program,  SMP,^*  was  developed  at  the  David  Taylor 
Naval  Ship  Research  and  Development  Center  (DTNSRDC)  in  1981  to  provide  a  standard 
ship  motion  prediction  tool  for  use  in  the  Navy's  ship  design  process.  This 
program,  also  known  as  SMP81,  provides  predictions  in  irregular  seas  of  the  six 
degree  of  freedom  ship  responses  (surge,  sway,  heave,  roll,  pitch,  and  yaw)  as  well 


*A  complete  listing  of  references  is  given  on  page  27. 
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as  of  the  longitudinal,  lateral,  and  vertical  responses  at  specified  locations  on 
the  ship.  Since  1981  a  number  of  corrections  and  improvements  were  made  to  SMP 
which  are  detailed  in  this  report. 

The  rodif ications  to  SMP,  known  collectively  as  198U  UPDATES,  include: 

1.  A  correction  to  bilge  keel  roll  damping  for  ships  that  have  both  a  bilge 
keel  and  e.  skeg  on  the  same  station. 

<?.  An  improvement  in  the  roll  prediction  using  a  different  theory  to  compute 
hull  lift  coefficients. 

3.  A  new  theory  to  compute  stabilized  ship  responses  using  active  antiroll 
fins  with  fixed  gains. 

U.  A  new  prediction  of  vertical  shear  force  and  vertical  bending  moment  at 
specified  cross  sections  on  the  ship. 

5.  The  lowest  frequency  of  encounter  computed  in  SMP  is  now  restricted  to 
0.05  radians  per  second  to  avoid  calculations  of  added  mass,  damping,  excitations, 
and  motions  at  near  zero  radians  per  second. 

6.  A  new  set  of  Severe  Motion  Tables  for  selected  responses  is  provided  as 
part  of  the  output. 

It  should  be  noted  that  no  major  alterations  were  made  to  SMP  in  either  input, 
output,  files,  or  the  manner  in  which  tile  program  is  executed.  The  updated  version 
of  SMP,  known  as  SMP8U ,  will  accept  input  deck3  prepared  for  SMPfil  as  well  as  files 
previously  generated  by  SMP81  using  these  decks.  There  are  changes,  however,  to 
both  the  input  deck  and  some  of  the  files  when  new  features  such  as  active  fins 
and/or  vertical  loads  are  selected. 

The  specific  changes  to  input,  output,  and  files  associated  with  active  fins, 
vertical  loads,  and  severe  net  ion  tables  are  described  ir.  separate  sections  of  this 
report . 

A  description  is  provided  of  the  method  used  to  rake  the  program  modifications 
to  SMP.  The  files  associated  with  these  program  changes  are  stored  on  a  disk  pack 
on  the  Control  Data  Corporation  CYBER  176  at  DTNSRDC.  The  last  section  of  this 
report  describes  these  files  and  details  the  procedure  used  to  retrieve  the  SMP84 
file.  A  listing  of  the  5MP0U  source  code  UPDATES  is  provided  in  Appendix  A. 

This  report  only  describes  the  changes  made  to  SMP01.  It  is  thus  meant  to  be 
used  as  a  supplement,  not  a  replacement,  for  the  SMP  User's  Manual.1  In  addition, 

p 

the  user  is  referred  to  Cox  and  Lloydc  for  more  details  concerning  the  design  basis 
for  using  active  anti roll  fins  in  roll  motion  stabilization. 
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SHIP  PARTICULARS 


Four  ships  are  used  for  illustrative  purposes,  these  are: 

1.  USCG  210-ft  Medium  Endurance  Cutter  (WMEC)  -  bilge  keel/skeg  damping 
correction 

2.  DE-10Q6  -  modification  to  hull  lift  terms 

3^  USCG  270-ft  Medium  Endurance  Cutter  (WMEC)  -  active  antiroll  fins 
4.  DDG-51  -  vertical  loads 

A  listing  of  the  hydrostatic  characteristics  of  each  ship  is  presented  in  Table  1 
and  Figure  1  contains  the  computer-drawn  underwater  hull  shapes. 


UPDATES  10  THEORY 

It  is  assumed  in  thi3  report  that  the  user  is  already  familiar  with  the  ship 
motion  theory,  variables,  coordinate  system,  files,  and  ir.put/output  schemes  that 
are  described  in  the  SMP  User's  Manual.^  These  details  will  not  be  repeated  here. 
Only  the  changes  made  to  SMP81  are  described  below. 


Dvr  r  /  r» i-r-w-i  MuinTlfn  ^Hnn/WT  .~\h1 

±u\J£m  ru^uL/ orvcA j  L/nrir  ahvj  uwiuuiV/ti.OM 


Most  of  the  problems  that  users  have  experienced  in  running  SMP  involve  the 
proper  selection  of  input.  Spline  routines  are  used  extensively  to  fit  the  hull. 
These  spline  fits  are  sensitive  to  curves  with  sharp  corners  or  that  have  uneven 
point  spacing  (some  points  too  close  together  and  other  points  spaced  too  far 
apart).  This  type  of  problem  is  corrected  by  using  a  different  point  selection  and 
thus  did  not  require  a  programming  change. 

A  different  problem  occurred  in  the  calculation  of  bilge  keel  viscous  roll 
Camping  for  ships  that  have  both  a  skew  and  a  bilge  keel  on  the  same  station.  A 
variable  called  the  "radius  of  the  bilge"  shown  in  Figure  2  was  incorrectly  com¬ 
puted  near  the  centerline  of  the  ship  instead  of  near  the  bilge  keel,  due  to  the 
presence  of  the  skeg.  This  error  caused  a  significant  increase  in  bilge  keel 
damping  even  though  the  skeg  may  have  only  overlapped  part  of  the  bilge  keel. 

This  error  has  been  corrected  in  SMP04.  An  exanple  illustrating,  the  change  in 
the  bilge  keel  component  of  roll  damping  is  shown  in  Figur>  3  for  a  210-foot  United 
States  Coast  Guard  cutter.  A  comparison  of  total  damping  between  the  SHPCl  and 
SMP84  predictions  as  well  as  measured  damping  from  full-scale  trials  is  also  shown 
in  Figure  3. 
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The  SMP&4  total  damping  preuiction  for  this  ship  reflects  an  additional 
decrease  in  roll  damping  due  to  a  change  in  the  method  of  computing  hull  lift  which 
is  discussed  in  the  next  section. 

MODIFICATIONS  TO  HULL  LIFT  TERMS 

One  of  the  improvements  made  to  SMP  in  1981  was  the  incorporation  of  hull 
and  appendage  lift  terms  in  the  lateral  equations  of  notion.  Validation  runs^ 
of  SMP  were  made  in  1981  for  the  DE-1006  which  showed  a  definite  improvement  in 
the  prediction  of  the  roll  transfer  functions  at  high  speed.  Although  the 
magnitudes  of  the  peak  values  were  predicted  quite  well,  there  wun  a  general 
tendency  to  shift  the  roll  transfer  functions  to  longer  wavelengths  than  showed 
by  experiment. 

Subsequent  checks  of  these  validation  runs  were  made  for  the  DE-lOOo  at  the 
lowest  GM  value,  GM3,  at  zero-speed  and  at  27  knots  in  beam  waves.  A  discrepancy 
was  found  between  the  zero-speed  undamped  natural  roll  period,  T^,  computed  in 
SMP  and  the  measured  roll  period  shown  in  Figure  t  (taken  from  Reference  3)  for 
the  GM3,  BKt  (no  bilge  keel)  condition. 

The  undamped  roll  period  was  computed  in  SMP  from  the  natural  roll  fre¬ 
quency,  tii,,  as 
v 

w^2  *  A  GM/(Ii^  +  A^)  (l) 


where  A  is  the  ship  displacement,  GM  is  the  metacentric  height,  1^  is  the 
structural  moment  of  inertia,  and  is  the  hydrodvnamic  moment  of  inertia. 
Here  Ij^  is  computed  as 


=  M  K, 


where  M  is  the  mass  of  the  ship  and  is  the  roll  radius  of  gyration. 

The  measured  value  of  provided  as  input  to  SMP  for  the  validation  runs 
was  35  percent  of  the  beam.  The  GM  value  computed  ir.  SMP  was  6  percent  of  the 


1* 


beam  which  agreed  vit.h  the  measured  value  of  GM.  The  roll  period  computed  in 

SMP  for  these  valusc  of  and  GM  was  13  seconds.  The  measured  roll  period 

shown  in  Figure  4  was  9*3  seconds.  This  discrepancy  in  period  was  resolved  by 

assuming  that  the  measured  value  of  K.  already  contained  c.  hydrodynamic  com- 

v 

ponent  as  well  as  a  structural  component,  i.e.,  the  natural  roll  frequency  should 
have  been  computed  for  this  ship  as 

w*  =  A  GM/IUU  (>*) 


Iu.  =  M  (K.  )d  (5) 

4t  ♦measured 

The  structural  radius  of  ©'ration,  required  as  input  for  3MP  was  then  recomputed  aa 


(K.  )2  =  (K  )2  -  AU/M  (6) 

which  gave  a  value  of  32  percent  of  the  beam.  The  nsw  roll  period  was  then  com¬ 
puted  by  SM°  as  9.3  seconds  which  agreed  with  the  measured  roll  period  at  zero 
knots. 

This  modification  to  the  radius  of  ©ration  corrected  paijt  of  the  shift  of 
the  predicted  roll  transfer  functions  to  longer  wavelengths  at  speed  but  did  not 
eliminate  it. 

Similar  roll  validation  runs  were  made  by  Schmitke1*  in  1978  for  the  DE-1006 
using  a  computer  program  called  ROLLRFT  which  also  incorporated  terms  for  hull 
and  appendage  lift.  These  validation  runs,  made  using  ROLLRFT,  snowed  a  similar 
good  agreement  in  predicting  the  peak  values  of  the  roll  transfer  functions  at 
high  speed.  In  addition,  the  roll  transfer  functions  were  predicted  quite  well 
at  each  frequency,  i.e,,  there  was  no  tendency  to  shift  the  transfer  functions 
to  longer  wavelengths  as  shown  in  SMP. 

A  comparison  was  made  of  the  theories  for  hull  and  appendage  lift  that  were 
implemented  in  SMP  and  ROLLRFT.  It  was  determined  that  the  only  significant 
difference  in  these  lift  theories  was  in  the  theory  used  to  obtain  hull  lift.  A 
discussion  of  the  differences  of  the  two  theories  of  hull  lift  as  implemented  in 
SMP81  and  ROLLRFT  is  provided  next. 
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The  hull  la  treated  In  SMP8l  as  a  low  aspect  ratio  lifting  surface  that 
generates  a  sway  lift  force  and  a  corresponding  roll  lift  moment  at  the  x,y,z 
location  of  the  center  of  pressure 


xcp  51  0 

(7a) 

ycp  = n 

(7b) 

zcp  -  '  (vcg  +  ^ 

(7c) 

where  V,_  is  the  vertical  center  of  gravity  of  the  ship  referenced  to  the 
waterplane  and  T  is  the  draft  of  the  ship. 

The  hull  lift  terms  that  are  aided  to  the  left-hand  side  of  the  lateral 
equations  of  motion  are 


b22L  = 
b2UL  * 
c26L  * 
b42L  “ 
bl*l*L  * 


Fz  sin  r 


-  F,_  y  sinr 

it 


-  V  b 


22L 


cu6L  -  7  b2UL 

c62L  *  V  b22L 
c6UL  *  7  b2UL 
b66L  *  (y2^2)  b"2L 


(8) 

(9) 

(10) 

(11) 

(12) 

/  1  *>\ 
\iJ  / 

(11*) 

(15) 

(16) 


where  the  angle  r  that  the  hull  makes  to  the  horizontal  axis  is  -90  degrees  „  y 
is  the  moment  arm  from  the  center  of  gravity  to  the  center  of  pressure,  V  is  the 
ship  speed,  and  wg  is  the  wave  frequency  of  encounter. 


— 
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The  hull  lift  corrections  to  the  sway  exciting  force,  roil  exciting  moment. 

* 

s 

and  yaw  exciting  moment  are  determined  at  the  center  of  pressure  as 

•  ' — 

JO 

»v 

F?l  “  f 2  9*nr 

(17) 

iX 

fi.l  -  f2  y 

(16) 

•  ■  :v 

m 

t." 

_ _ — 

••  i- 

F6l  -  (V/img)  F2L 

(19) 

*  j 

where 

... 

f2  3  Fz  “  sin  r  sin  u  exP^k  2cp) 

(20) 

Here  u  is  the  wave  frequency  and  k  is  the  wave  number.  The  unit  sway  lift  force. 

t 

F 

Z,  15 

Fz  *  (p/2)  A  V  Lcs 

(21) 

"f 

B 

where  p  is  the  mass  density  of  water,  A  is  the  area  (span  x  mean  chord] 

and  Lcs 

is  the  lift  curve  slope. 
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The  hull  is  also  treated  as  a  low  aspect  ratio  lifting  surface  in 

the  Schmitke 

ROLLRFT  program  which  uses  the  theory  of  Mandel^.  This  theory  assumes 

that  &  force 

■ 

"couple"  is  generated  by  the  hull  at  speed  with  no  net  sway  lift  force 

or  roll  lift 

moment.  A  yaw  moment  is  developed,  however,  due  to  the  force  couple. 

x_r  is  computed  at  the  centroid  of  area  of  the  hull  as 

.  V 

„  _  /nT n  dn 

(22) 

L.  , 

cp  JV* 

where  T_  is  the  sectional  draft.  y_n  and  z.,,  are  assumed  to  be  zero. 

The 

moment  arm  used  for  the  force  couple  is  computed  as 

r  • 

where  Cp  is  the  prismatic  coefficient  and  Lpp  i3  the  length  between 

L 

-  *  • 

<■  « 

perpendiculars. 

T 

r, 

►  _ 

- 

_  _ 

A 
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shown  in  Figure  5  for  the  DE-1006  in  beam  waves  at  0  and  27  knots  for  the  GM3,  BK1 
condition.  It  should  be  noted  that  the  SMP81  and  SMP3U  transfer  function  com¬ 
parisons  are  identical  at  0  Knots. 


ACTIVE  ANTIROLL  FINS 

Appendages  are  input  into  SMP  in  order  to  compute  their  contributions  to  roll 
damping  due  to  lift.  These  appendages  include  fins,  rudders,  sKegs,  and  propeller 
shaft  brackets.  The  appendages  used  in  SMP  are  considered  to  be  passive.  One  of 
the  improvements  made  to  SMP  was  the  incorporation  of  active  antiroll  fins  in  the 
lateral  equations  of  motion.  The  active  fins  use  fixed  gains  which  vary  with 
speed.  The  gains  as  well  as  controller  characteristics  are  provided  as  input  to 
SMP  by  the  user.  A  method  for  determining  these  gains  is  discussed  later  in  this 
section. 

The  force,  F,  and  moment,  M,  developed  by  the  active  motion  of  a  fin  are 
F  -  n  {-  VKCP  -  *a)8  ♦  13  y  -  f0V(xcp  -  x#)|B  ♦  f„V2B} 

=  nF*  (32) 


M  =»  r  x  F 
cp 


=  r__  x  nF* 
cp 


(33) 


where  F*  is  the  magnitude  of  the  force  due  to  the  active  fin  motion,  8.  The  dots 
used  in  Equation  (32)  stand  for  time  derivatives  so  that  8  and  B  are  fin  i-ate  and 
fin  acceleration  respectively. 

The  variable  V  in  Equation  (32)  is  the  ship  speed  and  xg  is  the  longitudinal 
location  of  the  fin  shaft  referenced  to  the  longitudinal  location  of  the  center  of 
gravity  of  the  ship.  The  longitudinal  location  of  the  center  of  pressure  of  the 
fin,  x  ,  is  assumed  to  be  at  the  fin  shaft  for  a  passive  fin  but  to  move  aft  of 
t.he  fin  shaft  for  an  active  fin.  The  longitudinal  distance,  xcp  ~  xs,  is  assumed 
to  be  l/6th  of  the  mean  chord  for  an  active  fin. 
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The  fir  added  mas  term,  A^,  la  defir.ed  as 

Ap  «  *p  s(c/2)2  (3U) 

where  p  is  the  mass  density  of  water,  s  is  the  span,  and  c*  is  the  mean  chord.  The 
term  fQ  is  defined  as 

r0  -  (p/2)  A  Lcs  05) 


where  A  is  the  planform  area  of  the  fin,  and  Lcs  is  the  fin  lift  curve  slope.  The 
position  vector,  ,  is  defined  as 


CP 


ix 


cp 


jy, 


cp 


♦  kz 


cp 


(36) 


where  xcp,  ycp,  and  zcp  are  the  locations  of  the  center  of  pressure  of  the  fin 
referenced  to  the  vertical  center  of  gravity  of  the  ship. 

The  vector,  it,  in  Equation  (32)  is  the  unit  vector  normal  to  the  fin,  defined 
as 


*■  +  *  ♦ 
n  *  10  -  J3inr  -*■  kcosf 


(3T) 


n 

i 


where  T  is  the  angle  that  the  passive  fin  makes  with  respect  to  the  horizontal 
axis.  The  unit  vector,  it,  the  angle  r,  and  the  position  vector  ?cp  are  shown  in 
Figure  u. 

The  fin  rate  term,  3ApVfi,  in  Equation  (32)  is  known  as  the  Magnus  Lift,  term.^ 
This  lift  term  occurs  because  an  active  fin  behaves  in  a  tanner  similar  to  a 
rotating  cylinder  in  a  flew.  The  general  form  of  the  Magnus  Lift  term  is 

FML  *  2,(Vt/V)(2rs)(pV2/2)  (38) 

where  the  tangential  velocity  at  the  trailing  edge  of  the  fin,  V^,  is 

Vt  ■  (3c7'*)0  (39) 
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V' 


”1 


The  radius  r  of  the  rotating  fin  cylinder  is 


r  »  c72 


Substituting  Vt  and  r  into  Equation  (33) 


2*(3c/l*)  0  c  a(p/2)V 

3l  wpsCc/2)2!  ve 


The  components  of  sway  force,  Fy,  roll  moment,  M^,  and  yaw  moment,  M^,  are 
obtained  by  substituting  the  expressions  for  ?cp  and  rf  into  Equations  (32)  and  (33) 

Fy  3  -  sinT  F*  (**2) 

M  *•  v  F*  ( **3) 


where 


M*  *  ~  Xcp  8inr  *** 


y  *  ycp  ♦  tcp  sinr 


Tne  fin  stabilized  lateral  equations  of  motion  are  obtained  by  adding  Fy,  H^, 
and  to  the  left-nand  sides  of  the  sway  equation,  roll  equation,  and  yaw 
emjat i on  *  respectively* 

The  commanded  mo  cion  of  the  fin,  3,  is  determined  using  a  control  law 
operating  on  the  roll  notion  of  the  ship  and  its  tine  derivatives  es 

0  -  Gy  [Xk/(XA  •  Xb)|  ♦  U6) 

where  Gy  i3  a  speed -dependent  gain  factor.  Xk,  XR,  and  Xb  are  complex  coefficients 
defined  an 


”  <*£2  K3)  i  <^2 
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Xa  *  ^al  “  ^;2  *3^  *  *£*2 


(48) 


(bi  - 


2  b,)  +  i 


where  K^,  Kg,  and  K3  are  fin  controller  coefficients  proportional  to  roll  accelera¬ 
tion,  roll  velocity,  and  roll  angle;  a^,  &2,  and  a^  are  fin  servo  coefficients;  and 
b^,  b2,  and  are  fin  controller  compensation  coefficients. 

The  coefficients  Gy  and  Kj ,  aj ,  and  bj ,  where  J«1  to  3,  are  required  an  input 
in  SMP.  Nominal  values  for  these  coefficients  were  obtained  from  Reference  2  and 
are  shown  in  Table  2. 

The  free-stream  fin  lift  curve  slope  in  Equation  (35)  is  computed  in  as 

Lcs  *  1.8  »  af  /  [cos  A(ae2  sec2  A+  U)1^  +  1,8]  per  raaian  ($C) 

where  is  the  effective  aspect  ratio  of  the  fin,  i.e.,  2s/c  and  A  la  the  sweep 
angle  of  the  fin  quarter  chord. 

The  user  can  optionally  input  an  effective  Hit  curve  slope  which  takes  into 
account  fin  performance  degradation  due  to  the  effects  of  hull  boundary  layer, 
fin /bilge  keel,  and  fin/fin  interference  effects.  The  user  is  referred  to 
Reference  2  for  methods  that  can  be  used  -o  determine  these  degradation  effects. 

One  degradation  effect  discussed  in  Reference  2  that  should  not  be  used  involves 

fin-induced  sway  and  yaw  motions.  Reference  2  used  this  particular  effect  tc 
account  for  sway  and  yaw  notions  in  a  one  degrec-of-f reedoaj  roll  equation  to  compute 
stabilized  roll  for  active  fins.  SMP  does  not  require  tills  degradation  errect 
because  it  uses  a  three  degree-of-f reedom  math  node],  to  compute  stabilized  sway, 
roll,  and  yaw  directly  using  the  fin  sway  force,  fin  roll  moment,  and  fin  yaw 
moment  from  K^uations  (42),  (43),  and  (4k). 

Finally,  if  the  user  does  not  have  specific  values  for  the  fin  gain  factors 
Gv,  it  jilx  be  necessary  to  find  theta  iteratively  by  making  a  nu^>er  of  runs  of  SrtP 

for  a  range  of  Gy  values  using  program  option  0PTN«3  and  0PTN*5  in  Data  Card  Set  2. 

The  user  is  referred  to  Reference  1  for  the  details  of  running  SMP. 

The  procedure  involves  applying  a  specific  fin  angle  limit  criteria,  based  or 
occurrence  of  fin  stock  strength,  stall,  cavitation,  and/or  reduced  noise  consider¬ 
ations  (see  Reference  2),  to  the  root  mean  square  fin  angles  output  by  SMP  at  the 
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heading  which  produces  the  worst  roll  motion  in  a  design  seaway.  It  should  be 
noted  that  the  significant  wave  height  for  this  design  seaway  should  be  input  by 
the  user  in  Data  Card  Set  14.  A  root  mean  square  statistic,  STATIS»1.0Q,  should 
also  be  input  in  Data  Card  Set  14  when  making  the  iterative  runs  of  SMP. 

Typically,  the  user  selects  a  fin  limit  angle,  B^IM’  (e*8*»  21  degrees  at  15 
knots)  for  each  speed  with  a  one  in  ten  probability  of  heing  exceeded.  This  means 
that  for  a  particular  speed,  the  root  mean  square  fin  angle,  Sjy^g,  is  determined  as 

Arms  ~  ^l) 

where  2.146  is  the  Rayleigh  constant  used  to  compute  this  probability  of 
exceedance. 

Next,  the  user  plots  the  fin  angle  RMS  values  at  the  worst  roll  heading  in  the 
design  seaway  as  a  function  of  the  G  values  that  were  used  in  the  various  runs  of 
SMP.  The  fin  RMS  limit  criteria,  for  each  speed  are  then  plotted.  The 

appropriate  value  of  Gy  for  each  speed  can  then  be  determined  from  the  intersection 
of  &t  that  speed  and  the  curve  of  computed  RMS  fin  angle  as  a  function  of  G. 

An  example  of  this  type  of  plot  is  shown  in  Figure  7  for  a  270-foot  United  States 
Coast  Guard  cutter, 

VERTICAL  LOADS 

The  theory  of  Salvesen,  Tuck,  and  Faltinsen1  is  use  to  compute  the  ver¬ 
tical  wave-induced  loads  at  specified  cross  sections  of  a  ship  advancing  at 
constant  forward  speed  with  arbitrary  heading  in  regular  sinusoidal  waver..  These 
loads  are  expressed  for  a  given  ship  speed,  heading  angle,  and  frequency  of 
encounter,  ut,  e.s  : 

Vj  =  tj  cos(u£t  +  5^  )  (52) 

where  bj  is  the  load  amplitude  with  J=3  referring  to  the  vertical  shear  force  and 
J=5  referring  to  the  vertical  bending  moment-  The  phase  angles  ij  refer  to  the 
phase  lead  of  load  with  respect  to  the  maximum  wave  elevation  at  the  origin  of 
the  x,y,z  coordinate  system  shown  in  Figure  8.  This  right-handed  coordinate 
system  is  moving  with  the  constant  mean  forward  speed  of  the  sh  j  vith  the  origin 
lying  in  the  undisturbed  free  surface  and  located  at  the  longitudinal  center  of 


gravity.  The  coordinate  system  is  defined  with  z  positive  vertically  upward 
through  the  center  of  gravity  of  the  ship,  y  positive  to  port, 'and  x  positive  in 
the  direction  of  forward  motion  of  the  ship. 

The  sign  convention  of  the  vertical  loads  is  also  shown  in  Figure  8.  The 
loads  are  located  at  the  specified  croso  section  with  the  vertical  shear  force, 

Vj,  positive  upward  and  the  vertical  bending  moment,  V^,  positive  bow  down.  It 
should  be  noted  that  is  actually  the  moment  about  the  horizontal  axia.  It  is 
referred  to  as  the  vertical  bending  moment  by  convention  since  it  is  the  moment 
due  to  the  vertical  forces. 

The  vertical  shear  force  is  computed  es  the  difference  between  the  inertia 
force  and  the  sum  of  the  external  forces  acting  on  the  portion  of  the  hull  forward 
of  the  specified  cross  section: 

v3  =  I3  -  (R3  +  Ej  ♦  d3)  (53) 

where  13  is  the  inertial  force,  R3  is  the  static  restoring,  is  the  exciting 
force,  and  D3  is  the  hydrodynamic  force  due  to  the  body  motion.  Similarly,  the 
vertical  bending  moment  is  computed  as  the  difference  of  the  BDment  due  to  the 
inertia  force  and  the  sum  of  the  moments  due  to  the  external  forces: 

V?  -  I5  -  (R5  +  Ej  +  d5)  (5**) 

The  vertical  inertia  force  is  equal  to  the  mass  times  the  acceleration: 

I3  *  /  n  (£3  -  nCj)  dn  (5.5) 

where  m  i3  the  section ,1  mas3  per  unit  length,  £3  is  the  heave  acceleration,  and 
is  the  pitch  acceleration.  The  sign  convention  of  the  six  degree  of  freedom 
motions,  ,  are  shown  in  Figure  8,  The  integration  is  over  the  portion  cf  the 
ship  forward  of  the  specified  cross  section. 

The  vertical  moment  of  inertia  is  defined  as: 

I5  -  -  /  »  («»-**)  (?3  -  n?5>  dn 


1*4 


(56) 


ft 

where  x  Is  the  longitudinal  location  of  the  cross  section  referenced  to  the  origin 
of  the  x,y,z  coordinate  system  in  Figure  8. 

The  vertical  hydrostatic  restoring  force  and  moment  over  the  portion  of  the 
ship  forvard  of  the  cross-section  at  x*  are  given  by: 

r3  *  -  <*  J  b  (43  -  nCj)  dn  (57) 

R5  “  06  /  h  (h»x*)  (£3  -  n^)  dn  (58) 

vhere  p  is  the  mss  density  of  water,  g  is  the  acceleration  of  gravity,  and  b  is 
the  sectional  beam. 

The  exciting  force  and  moment  over  the  portion  of  the  ship  forward  of  the 
cross  section  at  x*  are  defined  as: 

E3  «  pcA  {/  (f3  ♦  h3)dT,  ♦  |(V/ii^.)  h3)]frx*}  expUugt)  (59) 

E5  *  “  p;A  fKihx'l  (f,  ♦  h,)  ♦  (V/itt^)h3)dn  exp(ittpt)  (60) 

where  <;A  is  the  wave  amplitude  and  V  is  the  ship  speed. 

The  sectional  Freude-Krllofr  "force"  is  given  by: 

i 

f3  *  g  exp(-ikncos p)  Jq  N3  exp(ikysinp)  exp(kz)dl  (61) 

and  the  sectional  diffraction  force  is  given  by: 

h3  =  w  exp(-ikncos  p)  /c  (iN3  -  MgSinp)  exp(ikysinp)  exp(kZ)  ^  dl  (62) 

Here  k  is  the  wave  number  and  p  is  the  ship  heading  angle  relative  to  the  incident 
wave.  N2  and  N3  are  the  two-dimensional  sectional  normal  components.  Cn  denotes 
the  croBs-section  at  longitudinal  location 

The  hydrodynamic  force  and  moment  due  to  the  body  motion  on  the  portion  of  the 
ship  forward  of  the  cross  section  at  x*  are  given  by: 
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^3  *  “  /  {&33^3  ~  +  b33^3  "  **£5) 

"  (V/«ig^)b32^  +  ',a33^5* 

~  [Va33  (£3  -  ^$5)  ~  W«g2)b33(£3  -  nSj) 

-  (V2/we2)(r33^5  +  b33C5)]rT.x*  (63) 

and, 

D5  ~  /  ( n-x  ^  (a33^3  “  b^'5)  +  b33(  C3  -  nt^)  }  drt 

+  /  {Va33 (  C3  -  x*£5)  -  iv^h  ‘  x*%) 

“  ^^<‘£^^a33  £5  +  t>33 €5 )  }  dn  (6U) 

where  a33  and  b33  are  the  two-dimensional  sectional,  added-m&ss  and  damping  for 
heave . 

RESTRICTION  ON  ENCOUNTER  FREQUENCY  CALCULATION 

The  frequency  of  encounter,  wg,  is  computed  as 

«E  ■  |«-  (w2V/g)cosu|  (65) 

where  V  is  the  mean  forward  speed  of  the  ship,  u  is  the  heading  angle,  and  «  is  the 
wave  frequency. 

The  two-dimensional  velocity  potentials  as  well  as  the  sero  speed  added  mass 
and  damping  coefficients  are  calculated  in  SMP  over  a  fixed  rmngc  of  10  encounter 
frequencies  from  0.05  to  10.0  radians  per  second.  The  2-D  velocity  potentials  and 
added  masa  and  damping  coefficients  are  then  spline  fitted  over  thin  range  of  10 
encounter  frequencies.  The  2-D  velocity  potentials  and  added  Bass  and  damping  for 
specific  wg  values,  computed  in  Equation  (65)  fo;  various  ship  speeds,  headings  and 
wave  frequencies,  are  then  obtained  by  interpolation  from  the  spline  fitted  poten¬ 
tials  and  zero  speed  added  mass  and  damping  coefficients. 

The  spline  interpolation  routines  in  SMP  do  not  allow  extrapolation  outside 
the  range  of  the  independent  variable,  i.e.„  the  range  of  encounter  frequencies 
from  0.05  to  10.0  radians  per  second.  However,  it  was  found  that  at  high  speed  in 
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quartering/following  waves  that  values  near  zero  radians  per  second  were  com¬ 
puted  using  Equation  (65)  and  that  the  2-D  potentials  and  added  nass  and  damping 
coefficients  were  being  extrapolated  to  these  near  zero  encounter  frequencies.  The 
extrapolation  was  done  because  other  source  coding  was  used  for  interpolation 
instead  of  the  spline  interpolation  routines.  This  resulted  in  very  large 
unrealistic  values  for  the  wave  excitations  on  the  light  hand  side  in  the  equations 
of  motion  for  these  near  zero  encounter  frequencies. 

The  wg  values  computed  using  Equation  (65)  in  SMPB4  are  set  to  0.05  if  they 
are  less  than  0.05  radians  per  second.  In  addition,  a  separate  range  of  100 
encounter  frequencies  used  in  SMP  to  interpolate  encounter  spectra  was  modified  so 
that  the  lowest  encounter  frequency  would  be  0.05.  The  modified  set  of  encounter 
frequencies  is  now  0.05  (O.Oi)  O.58,  0.60  (0.02)  1,  1.1  (0.1)  2,  2.2  (0.2)  4,  4.4 
(0.4)  6.  The  variable  frequency  increment  is  shown  in  the  parenthesis. 

SMP  PROGRAM  CHANGES 

INPUT 

The  inp^*  "or  SMP  consists  of  hull  form  data,  loading  data,  appendage  data, 
point  location  data,  and  environment  data.  This  input  is  broken  down  into  15  Data 
Card  Sets  which  are  described  in  Appendix  C  of  Reference  1.  The  modifications  to 
SMP  to  incorporate  vertical  loads  required  changes  to  Data  Card  Set  2  (Program 
Options)  and  Data  Card  Set  6  (Undervater  Hull  Geometry).  Modifications  for  active 
ar'.iroll  fins  required  changes  to  Data  Card  Set  11  (Fin).  The  changes  to  these 
three  data  card  sets  are  described  below; 

Data  Card  Set  2,  Program  Options 

A  new  printing  option  (LRA0PR,  integer,  column  25)  is  provided  to  print  out 
the  vertical  shear  force  and  vertical  bending  moment  response  amplitude  operators 
(RAO)  and  phase  angles.  A  new  load  RAO  file  (TAPL10)  is  generated  only  when  the 
DRA0PR  option  is  selected.  The  values  of  LRA0PR  are 

0  or  blank  -  No  load  RAO  printout. 

1  ~  Print  out  the  load  RAO's  and  generate  a  load  RAO  file  when  0PTN=2 
through  6.  0PTN  is  the  rtfijor  program  option  specified  in  column  5 
of  this  data  card  set.  A  loac:  coefficient  file  (TAPEU)  must  be 
attached  when  0PTN=4  through  6.  An  origin  file  (TAPEll)  must  also 
be  attached  when  0PTN=6. 


In  addition,  the  load  variable  NLOADS ,  mat  be  selected  in  Data  Card  Set  6  to 
obtain  the  load  RAO  printout. 

Data  Card  Set  6,  Underwater  Hull  Geometry 

A  new  option  (NLOADS,  integer,  columns  9-10)  is  provided  on  the  first  card  of 
this  data  card  set.  NLOADS  specifies  the  number  of  stations  (maximum  of  10)  where 
vertical  loads  eu-e  to  De  computed.  No  vertices,  loads  are  computed  when  NLOADS  is 

specified  as  either  0  or  blank.  If  the  load  option  is  selected  (NLOADS  >0),  two 

additional  sets  of  information  mist  be  provided  at  the  end  of  this  data  card  set 

after  the  stations  and  hull  offsets  or  Lewis  forms  have  been  input. 

First,  the  static-  weight  ( SWGHT ,  real  array,  8F10.4)  is  input  for  each  of  the 
stations  specified  previously  in  this  data  card  set.  The  weight  units  are  metric 
tons  (mass  unit)  if  PUNITS*"MET1£R"  or  long  tons  (weight  unit)  if  HJINITS*"FEETh. 

The  variable  PUNITS  is  specified  in  Data  Card  Set  3.  The  weight  curve  for  the 
DDG-51  is  shown  in  Figu: .  9. 

Second,  the  stations  (X.LDSTN,  real  array,  8F10.4)  where  loads  are  to  be  calcu¬ 
lated  arc  input.  The  variable  NLOADS  determines  the  number  of  load  stations  that 
are  specified.  A  load  station  must  correspond  exactly  to  one  of  the  station  num¬ 
bers  specified  earlier  in  this  data  card  set. 

Data  Card  Set  11,  Fin 

Two  nav  variables  associated  with  active  fins  are  provided  on  the  first  card 
of  this  data  card  set.  The  first  variable  (IACTFN,  integer,  column  10)  specifies 
whether  the  fins  are  active  (IACTFN=1)  or  passive  (IACTFN^O  ).  The  second  variable 
( TFCLCS .  integer,  column  15)  allows  the  user  (IFCLCS*1)  to  input  an  effective  lift 
curve  slope  for  each  3peed  and  fin. 

Four  new  cards,  denoted  as  cards  1.1,  1.2,  1.3,  and  1.4,  oust  be  input  after 
the  first  card  when  active  fins  ewe  selected.  These  cards  are  skipped  when 
:actfn=o. 

Card  1,1  -  FORMAT  (8F10.4) 

(1)  FGAIN  (array),  real,  columns  1-10,  11-20 . l(NVK-l)»10+l|  -  lNVK*10l  , 

speed-dependent  fin  gain  factors,  Gy.  NVK  is  the  number  of  ship 
speeds. 
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Card  1.2  -  FORMAT  (3F10.4) 

(1)  FK  (array),  real,  colunms  1-10,  11-20,  21-30,  Tin  controller  coefficients 
where  PTC(l)  i3  proportional  to  roll  angle,  re(2)  io  proportional  to  roll 
velocity,  and  FK(3)  ia  proportional  to  roll  acceleration. 

Card  1.3  -  FORMAT  (3F10.4) 

(l)  FA  (array)  real,  columns  1-10,  11-20,  21-30,  fin  servo  coefficients 
Card  1,4  -  FORMAT  (3F1G.4) 

(1)  FB  (array),  real,  columns  1-10,  11-20,  and  21-30,  fin  controller 
compensation  coefficients. 

Nominal  values  for  FK,  FA,  and  FB,  taken  from  Reference  2,  are  shown  in 
Table  2.  FGAIN  is  either  known  for  existing  ships  or  determined  by  making  itera¬ 
tive  runs  of  SMP  following  the  procedure  described  in  the  active  fin  section  of 
this  reqport . 

The  next  card,  denoted  as  1.5,  is  required  for  each  fin  set  if  the  user  wants 
to  input  (IFCLCS=l)  and  affective  lift  curve  slope.  Card  1.5  follows  Card  1.4  when 
1ACTFN=1  or  Card  1  when  IACTFN=0.  Card  1.5  Is  skipped  when  IFCLCS=«Q. 

Card  1.5  -  FORMAT  (8F10.4) 

(l)  FCLCS,  array,  real,  columns  1-10,  11-20,  .  .  1  ( NVJC-1  )e10+lJ  -  (NVKe10]  , 

speed-dependent  effective  fin  lift  curve  slope,  Lca,  for  a  particular 
fin  set.  The  user  is  referred  to  Reference  2  as  well  a a  the  section 
on  active  fins  for  methods  than  can  be  used  to  determine  LCB. 

Note  that  Card  1.5  is  input  consecutively  for  each  fin  set  prior  to  inputting  Card 
2,  which  provides  the  geometric  description  of  the  fins. 

OUTPUT 

The  basic  format  of  the  SMP  output  remains  unchanged  in  SMP84.  However,  the 
following  sections  were  modified  to  provide  output  for  active  antiroll  fins  and/or 
vertical  loads: 

1.  Input  Card  Description  for  fins  and  loads. 

2.  Response  Amplitude  Operators  for  loads. 

3.  Response  Statistical  Value  (RSV)/Encountered  Modal  period  (Tqjj.)  tables 
for  fins  and  loads. 

Two  examples  ore  provided  to  illustrate  these  output  codifications.  The  first 
example  is  for  active  antiroll  fins  on  a  270-foot  United  States  Coast  Guard  Cutter. 
The  SMP  input  deck  for  this  ship  is  shown  in  Table  3.  Table  4  shows  the  change 
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in  the  input  card  de8cription  for  Data  Card  Set  11  to  print  out  the  values  of  the 
variables  IACTFN,  IPCLCS,  FGAIN ,  FK,  FA,  FB.  and  FHCLCS.  An  effective  lift  curve 
slope  is  provided  in  the  speed-dependent  input  variable  FHCLCS  to  aULov  for  the 
effect  of  hull  boundary  layer  aa  veil  as  the  fin  performance  degradation  due  to  the 
presence  of  a  bilge  keel  aft  of  the  fin.  The  procedure  used  to  compute  the  effec¬ 
tive  lift  curve  slope,  taken  from  Reference  2,  is 


(Lcs^E  *  I'bl^1  "  *bk^Lc8 

where  Lcs  is  the  free  stream  fin  lift  curve  slope  defined  in  Equation  (50)»  The 
hull  boundary  layer  correction,  e^,  is  defined  as 

ebi  »  1  -  0.50(«/s)/Lcb  (67) 

where  s  is  the  fin  span.  The  boundary  layer  thickness,  4,  is  defined  as 

4  =*  0.377  xFp(Re)-°*2  168) 


where  xFp  is  the  distance  between  the  ship  forward  perpendicular  and  the  longitudi¬ 
nal  location  of  the  center  of  pressure  of  the  fin,  (xCp)fin»  The  Renolda  number. 


Re,  is 


*  XppV/ v 


(69) 


where  V  i3  the  ship  speed  and  v  is  the  kinematic  viscosity. 

The  bilge  keel/fin  degradation  correction,  e^,  is  defined  aa 

ebk  »  0.22  (ae)bk  (l  ♦  il  +  (s/d^)2!^2}  (Abk/s2)  (7u) 

where  A^  is  the  area  of  the  bilge  keel  and  d^b  la  the  distance  between  (*Cp)fin  and 
the  midlength  location  of  the  bilge  keel,  (xCp)bk»  The  aspect  ratio  of  the  bilge 
keel,  (ae)bk»  is  defined  as 

(ae>bk  *  2  sbk/cbk  (71) 


where  and  cbk  are  the  bilge  keel  span  and  mean  chord  respectively. 
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The  kinematic  viscosity,  y,  is  provided  as  part  of  the  SMP  input  in  Data  Card 
Set  2.  The  values  for  Abk,  3bk,  c^,  s,  Lc8,  (xcp)fin  and  (xcp)bk  are  printed  out 
in  the  Roll  Damping  Tables  which  can  be  obtained  by  specifying  RLDMPR=1  in  Data 
Card  Set  2.  The  reference  for  (xcp)f.in  and  (xCp)fck  is  the  longitudinal  center  of 
gravity  of  the  ship,  LCG.  LCG  itself  is  referenced  to  the  forward  perpendicular 
and  is  printed  out  in  the  Hydrostatic  Table. 

The  values  computed  for  e^,  and  (Lca)E  for  the  USCG  2T0-foot  cutter  were 
0.9684,  0.1552,  and  2.118,  respectively,  at  a  ship  speed  of  15  knots.  Thus  the 
value  used  for  the  input  variable  FNCLCS  at  15  knots  was  2.118.  The  uncorrected 
free  stream  value  of  2.529  was  used  for  FNCLCS  at  zero  knots. 

It  is  necessary  for  the  user  to  make  an  initial  run  of  SMP  for  the  unstabi¬ 
lized  ship  (fin  gain  of  zero)  using  the  no...  motion  only  option  in  Data  Card  Set  2 
with  RLDMFR=1  in  order  to  obtain  the  Hydrostatic  Table  and  the  Roll  Damping  Table. 
The  values  of  FNCLCS  can  then  be  computed  fro:  tnfcrxation  provided  ir.  these  tables 
and  provided  in  the  SMP  input  for  succeeding  3u2  runs  for  active  fin3. 

The  values  used  for  FK,  FA,  and  FB  are  taken  from  Reference  2.  The  values  for 
PGAIN  were  obtained  using  the  iteration  procedure  described  in  the  section  of  this 
report  on  active  fins.  A  fin  limit  angle  of  21  degrees  was  applied  at  a  ship  speed 
of  15  knot3,  heading  of  105  degrees  (head  seas  equals  0  degrees),  in  short-crested 
seas  with  significant  wave  height  of  3*96  meters  ( 13  feet)  anft  modal  wave  period  of 
9  seconds.  The  nondimens ional  roll  decay  coefficients  are  shown  in  Table  5.  The 
short-crested  RSV/TQE  printout  for  unstabilized  roll  single,  stabilized  roll  angle, 
stabilized  roll  velocity,  fin  angle,  and  fin  velocity  are  shown  in  Tables  6  through 
10.  Note  that  the  significant  single  amplitude  statistic  (STATISTS,  Data  Card  Set 
14)  was  used  in  thi3  example. 

The  second  example  illustrates  the  vertical  load  output  for  the  DDG-51.  The 
SMP  input  deck  for  this  ship  is  shown  in  Table  11,  Table  12  shows  the  change  in 
the  input  card  description  for  Data  Card  Set  2  to  print  out  the  value  for  LRAOPR. 
LRA0PR=1  in  this  example  specifies  that  the  Response  Amplitude  Operators  for 
Vertical  Shear  Force  and  Vertical  Bending  Moment  are  to  be  printed  out.  Table  13 
shows  the  change  in  the  input  card  description  for  Data  Card  Set  6  for  loads.  The 
variable  NLOADS  on  the  first  card  of  this  card  set  specifies  that  the  vertical 
loads  for  one  station  are  to  be  computed.  At  the  end  of  this  card  set,  the  weight 
curve  in  units  of  long  tons  is  printed  out  as  a  function  of  station.  The  specific 
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station,  10.5  in  this  example,  vmere  vertical  loads  are  to  be  calculated  is  printed 
out  next  after  the  weight  ci rve. 

Table  lU  shows  the  RAO's  *;.d  phase  angles  for  the  Vertical  Shear  Force, 
called  V.SHEARlVj)  on  the  p.  Int.-i’t,  and  Vertical  Bending  Moment,  called  V.M0M.(V5), 
at  a  ship  speed  of  is  knots  and  a  heading  angle  of  L5  degrees.  This  Is  the  same 
format  used  to  print  out  the  RAO's  and  phases  for  the  six  degree  of  freedom  notions 
in  SMP. 

Table  15  shows  the  RSV/Tq^  values  for  the  Vertical  Shear  Force,  called 
V. SHEAR  FORCE  in  this  printout,  at  station  10.5  in  short-crested  seas  for  a  signi¬ 
ficant  wave  height  of  10  feet.  The  force  physical  units  are  in  long  tons/100 . 

Table  16  shows  the  RSV/TqE  values  for  the  Vertical  Bending  Moment,  called  V.BEND. 
MOMENT,  for  the  same  sea  condition.  The  moment  physical  units  are  in  foot-long 
tons/ 10000.  The  format  for  these  load  RSV/TqE  values  is  Identical  to  that  used  for 
the  other  ship  responses  printed  out  by  SMP. 

A  new  output  section  which  provides  tables  for  severe  notions  is  discussed  in 
the  next  section  of  this  report. 

SEVERE  MOTION  TABLES 

For  design  purposes  it  is  important  to  know  the  worst  (maximum)  values  of  the 
most  important  ship  responses*  heave,  pitch,  sway,  roll,  yaw,  £nd  the  vertical  and 
lateral  accelerations  for  up  to  four  locations  on  the  ship  specified  by  the  user. 
The  new  severe  motion  tables  provide  this  information  in  both  long -crested  and 
short-crested  seas  for  up  to  four  seaways.  Each  seaway  is  defined  by  a  significant 
wave  height  input  by  the  user  and  a  most  probable  wave  period  determined  in  SMP. 
Table  17  shows  the  most  probable  periods  used  in  SMP  for  various  ranges  of  signi¬ 
ficant  wave  height  given  in  meters. 

Each  severe  motion  table  is  organised  into  two  parts.  The  first  part  provides 
the  maximum  Response  Statistical  Value  (RSV)  and  associated  Encountered  Modal 
period  (TqE)  for  each  of  the  responses  listed  above.  The  ship  speed  and  beading 
where  these  maximum  responses  occur  are  also  given  in  this  first  part.  The  maximum 
RSV/TqE  values  are  obtained  from  the  standard  response  statistical  tables  that  are 
output  by  SMP.  It  should  be  noted  that  the  statistic  used  for  these  tables  (RMS, 
SIGNIFICANT  SINGLE  AMPLITUDE,  etc.)  i3  input  by  the  user. 

Another  set  of  information  that  is  useful  to  know  is  what  the  associated 
responses  are  at  the  conditions  (speed  and  heading)  where  the  mximum  responses 
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occur.  This  information  is  provided  in  matrix  form  in  the  second  part  of  each 
severe  motion  table.  Each  rov  of  the  matrix  is  associated  with'  one  of  the  maximum 
responses  which  is  listed  at  the  beginning  of  the  row.  The  speed  and  heading  where 
the  naximum  response  occurs  is  listed  next  in  the  row.  The  ISV-Tqj.  values  for  all 
of  the  responses  at  this  speed  and  heading  are  listed  next  in  the  row. 

An  example  of  the  new  severe  motion  table  is  provided  in  Table  18  for  a 
270-foot  United  States  Coast  Guard  Cutter. 

FILES 

Three  new  files  have  been  added  to  the  SMP  and  one  existing  file  has  been 
modified.  The  new  new  files  are:  (l)  Load  Coefficient  file  (LCOFIL.TaPEU) ;  (2) 
Load  Response  Amplitude  Operator  file  ( LRAFIL ,TAPE10 ) ;  and  (3)  Severe  Motion  file 
(SEVFIL.TAPEll*).  The  existing  file  that  was  changed  is  the  FMS/Tq£  file 
( RMSFIL .TAPE13 ) •  A  descripton  of  each  file  as  to  its  contents  and  where  it  is 
generated  and/or  accessed  in  SMP  is  presented  next.  It  should  be  noted  that  files 
in  SMP  are  identified  by  both  name  and  number,  me  user  can  catalog  and/or  attach 
files  only  by  using  the  tape  numbers, 

(1)  Load  Coefficient  file  -  LCOFIL  (TAPE!*) 

This  file  is  generated  in  subroutine  OOFOUT  when  the  program  option 
(OPTN  in  Data  Card  Set  2)  is  either  2  or  3  and  the  load  option 
( NLOAD  in  Data  Card  Set  6)  is  made  greater  than  zero,  mis  fixe 
must  be  attached  wen  the  user  selects  OPTN  greater  than  3  and 
NLOAD  is  greater  than  zero. 

The  file  contains  the  sectional  heave  exiting  force,  added-mass, 
and  damping  for  each  station  on  the  ship,  mis  file  is  read  in 
subroutine  RAOPHS  for  the  rms/toe  calculations  and  in  subroutine 
LRAOUT  to  print  out  the  load  response  amplitude  operators. 

(2)  Load  Response  Amplitude  Operator  file  -  LRAFIL  (TAPE10) 

This  file  Is  generated  in  subroutine  LRAOUT  and  contains  the  response 
amplitude  operators  for  the  vertical  shear  force  and  vertical  bending 
moment  for  load  stations  selected  by  the  user.  This  file  can  be 
cataloged  by  the  user  for  the  purpose  of  transferring  this  information 
to  some  other  computer  program. 

(3)  Severe  Motion  file  -  SEVFIL  (TAPElU) 

This  is  a  random  access  file  generated  in  subroutine  RMSOUT  that 
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contains  the  Response  Statistical  Values/Encountered  Modal  Periods 
(RSV/TQg)  for  heave,  pitch,  sway ,  roll,  yaw,  and  the  vertical  and 
lateral  accelerations  for  up  to  four  point  locations  specified  by  the 
user.  This  file  is  read  in  subroutine  SEVMOT,  which  prints  out  the 
Severe  Motion  Tables, 

(4)  RMS/TQE  file  -  FMSFIL  (TAPE13) 

This  file  contains  the  root  mean  square  (RMS)  values  and  the  encountred 
modal  periods  for  all  motions  defined  by  the  user  in  the  input  to 
SMP.  This  file  already  exists  in  SMP  but  was  modified  to  add  new 
responses  for  fins,  fin  angle  and  fin  velocity,  as  well  aa  the  vertical 
shear  force  and  vertical  bending  moment  for  up  to  10  stations  specified 
by  the  user. 


SOURCE  CODE  MODIFICATIONS 


The  Standard  Ship  Motion  Program,  SMP,  was  written  in  FORTRAN  IV  for  the 
CONTROL  DATA  CORPORATION,  CDC,  computers  at  DTNSRDC.  The  listing  cf  the  SMP  source 
code  was  provided  in  Appendix  I  of  Reference  1.  Each  line  in  this  listing  contains 
80  columns.  The  FORTRAN  source  code  is  contained  in  the  first  72  columns.  In 
addition,  a  subroutine  name  and  sequence  number  are  provided  in  ccluwis  73  through 

8o. 

A  batch  editor  called  UPDATE®  was  used  to  assemble  the  FORTRAN  source  code, 
subroutine  names ,  and  sequence  numbers.  The  same  editor  was  used  to  generate  the 
changes  to  SMP  described  in  this  report.  A  listing  of  these  changes  is  provided  in 
Appendix  A. 
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code.  The  editing  commands  are  Identified  by  an  asterisk  in  column  1.  The  editing 
commands  reference  specific  lines  in  the  SMP  listing  by  the  subroutine  names  and 
sequence  numbers  associated  with  these  lines.  A  description  of  these  update 
editing  comrands  follows: 

•ID  SMP84  -  The  name  SMP84  is  given  to  the  set  of 

modifications 

*1  "NAME" . N  -  The  FORTRAN  Instructions  which  follow  this  command 

are  inserted  in  subroutine  "NAME"  after  line  N 

*D  "NAME" .N( ,M)  -  Line  N  (optionally  lines  N  through  M)  in 

subroutine  "NAME"  are  to  be  deleted.  Any 
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FORTRAN  instructions  which  follow  this 
command  are  inserted  after  the  deletion 

*AF  -  New  FORTRAN  subroutinet  -nich  follow  this  command 

are  appended  to  the  end  of  the  SMP  source  code 

*DECK  "NAME"  -  Name  to  be  assigned  to  the  Fortran  instructions 

which  follow  tliis  card 

A  listing  of  the  segmentation  cards  required  to  load  SMP84  is  provided  in 
Table  19. 

DISK  PACK  STORAGE  OF  SMP  UPDATES 

Permanent  storage  of  SMP  and  its  UPDATES  is  maintained  at  DTNSRDC  on  the  CDC 
CYBER  1?6  computer.  Access  to  SMP  and  related  files  can  be  obtained  by  attaching  a 
disk  pack  and  copying  the  desired  file  to  the  main  disk  memory.  The  following  is 
an  example  of  the  CDC  control  cards  required  to  retrieve  SMP84  from  the  disk  pack: 

CHZM.CM55000 ,T5,P3. 

CHARGE  ,CHZM , XXXXOOOOO . 

PAUSE.  JOB  REQUIRES  DISK  PACK  DVUjOl 

MOUNT,  VSN*DVl*901  ,SN=TAPK06. 

ATTACH , A , SM  P0  4  ABSOLUTE , ID-CHZM ,SN»TAPK06 . 

REQUEST, SMP84,*PF. 

C0PYE,A,SMP84. 

CATALOG ,SMP84  ,SMP84ABSCLUTE,ID*=CHZM,MR=1 . 

EOF 

Examples  of  CDC  control  cards  required  to  run  various  options  of  SMP  can  be 
found  in  Appendix  E  of  Reference  1.  The  core  requirement  to  run  SMP 8 4  is  now 
150000  octal.  This  represents  a  change  from  100000  octal  required  to  run  SMPOl ,> 

The  CDC-6700  computer  at  DTNSRDC  has  been  replaced  by  &  CDC  CYBER-176  computer 
which  is  10  times  faster  than  the  CDC-6700.  Thus  the  SMP  run  times  specified  in 
Appendix  E  of  Reference  1  should  be  reduced  by  a  factor  of  10. 

The  following  files  associated  with  the  SMP  UPDATES  are  stored  on  the  disk 
paca  and  are  available  to  the  user: 

1.  SMP84ABS0LUTE  -  The  absolute  version  of  SMP84  used  for  production  runr/ng. 

2.  SMP04  CHANGES  -  The  file  containing  the  1984  UPDATE  modifications  to 

SMP81. 
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3.  SMP8 1* UPDATE  -  The  UPDATE  file  containing  the  original  SMP81  source 

code  as  veil  as  the  I98U  UPDATES, 

1*.  SMP6UC0KPILE  -  The  source  file  for  SMP8!». 

5.  SMP8U0BJECT  -  The  object  file  for  SMP8L. 

6.  SMP84SEGCABDS  -  The  segmentation  loader  cards  for  SMF0L. 

T.  SMP8I1LIBRARY  -  The  library  of  object  subroutines  for  SMP84. 

These  seven  files  are  stored  on  Disk  Pack  TAPK06  and  user  ID*CHZM.  The  set- 
name  (SN)  for  this  disk  pack  is  TAPK06  and  the  Volume  Serial  Number  (VSN)  is 
DVUcjOl.  All  seven  files  can  be  copied  to  the  main  disk,  memory  using  the  COPYE 
utility. 

For  further  information  regarding  CX  control  cards,  disk  packs,  and  copy 
utilities,  the  user  is  referred  to  the  DTNSRDC  Comjuter  Center  CDC  Reference 
Manual^ . 
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Figure  2  -  Example  of  the  Determination  of  the  Radius  of  the  Bilge 
Circle  in  SMP81  and  SMPdL  for  a  Station  where  the  Bilge 
Keel  Overlaps  the  Skeg 
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Figure  3a  -  Bilge  Keel  Component  of  Roll  Damping 


a;  Figure  3b  -  Total  Roll  Damping 


Figure  3  -  Example  of  the  Change  in  Roll  Damping  Coefficient  Due  to 
*-  the  Bilge  Keel/Skeg  Correction  for  the  USCG  210-ft  Medium 

Endurance  Cutter  (WMEC)  for  a  Ship  Speed  of  16  Knots 
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Figure  5  -  Coaparison  of  Predicted  (3MF84  and  SMPdl)  and  Meaeured 
Hor.dimensional  Roll  Transfer  Functions  for  the  DE-1006  at 
0  and  27  Knots  :Ln  Beam  Waves  for  GM3,  Bid 
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Figure  8a  -  Sign  Convention  of  Trans latory  and  Angular  Displacements 
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Figure  8b  -  Sign  Convention  of  Dynamic  Wave-Lead  Components 


Figure  8  -  Sign  Conventions  of  Motions  and  Loads 
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TABLE  1  -  SHIP  PARTICULARS 


USCG 

210-ft. 

WMEC 

DE-1006 

USCG 

270  ft 
»CC 

DD0-51 

Length  (Lpp),  meters 

6l.O 

93.9 

77. 7 

11*2.0 

Beam,  (E),  meters 

1C.1 

11.0 

11.6 

18.0 

Draft  (T),  meters 

3«3 

3.7 

U.2 

6.3 

Displacement  ( A)  S.W. , 
metric  tons 

1025.2 

1Q5U.2 

1818.7 

8560.2 

Mecacentric  Height  (GM),  %B 

6.1 

6.3 

8.2 

8.1 

Center  of  Gravity  (KG) , 

1*5.5 

1*3.5 

1*1*. 5 

1*1.3 

LCG*,  *LPP 

51.6 

51. T 

51.3 

50.1 

Natural  Roll  Period,  (T$), 
seconds 

10.? 

9.3 

IO.U 

13.1* 

Pitch  Gyradius  (Kg),  >Lpp 

25.0 

25.0 

2h.O 

25.0 

Roll  Gyradius,  (K$),  %B 

38.0 

32.1 

1*1.1 

1*0.0 

Yaw  Gyradius  (K^),  %Lpp 

25.0 

25.0 

25.0 

25.0 

Block  Coefficient  (CBj 

0.1*9 

0.51 

0.1*T 

0.52 

Bilge  Keel  Length,  <Lpp 

32.5 

29.6 

20.1* 

20.3 

Bilge  Keel  Span,  meters 

0.59 

0.1*7 

0.61 

0.9! 

Fin  Area,  j,.  meters 

— 

2.33 

.m 

•Referenced  to  F.P 


TABLE  2  -  NOMINAL  VALUES  OF  FIN  CONTROLLER  COEFFICIENTS,  K  j .  FIN  SERVO 
COEFFICIENTS,  aj ,  AND  FIN  CONTROLLER  COMPENSATION  COEFFICIENTS,  bj 


Kl 

»  1,0 

a^  *  1.0 

*1 

■  J.C 

k2 

*  2.5 

a2  91  Oj16 

b2 

«  0.025 

K3 

-  1.0 

aj  ■  0.63 

*3 

•  0.0y2 
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SMP  INPUT  DECK  FOR  USCG  270-FT  MEDIUM  ENDURANCE  CUTTER  (WMEC) 
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TABLE  13  -  EXAMPLE  OP  INPUT  CARD  DESCRIPTION  FOR  DATA  CARD  SET  6  FOR  LOADS 
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TABLE  17  - 

DETEKMIHATIOH  OF 

MODAL  WAVE  PERIODS  AS  /  PU1CTI0I  OF 
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SIGSIFICAHT  WAVE  HEIGHT 

FOR  THE  SEVERE  MOTIOB  TABLES 
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7-21  seconds  for 


the  range  of  sod el  ware  periods  used  in  SMP: 
ships  with  roll  periods  <_  15  seconds 
ships  with  roll  periods  ~  15  seconds 
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TABLE  19  -  LISTING  07  SMP84  SEGMENTATION  CARDS 


TREE  SMPR1-(  INPUT- (READ, HSTAT,  HSTOUT)  ,  REGVAV- (HYOCAL- (HY030 ,  T30AI© . 
. COFOUT ) . RQSASE . EQMOTN) . I  ROSE * . OUTPUT  - ( RAOOUT . LRAOUT . RIISOUT ) ) 

SNPSi  INCLUDE  SNP81 
INPUT  INCLUDE  INPUT 

READ  INCLUDE  READ , AINPUT , GENOFS . BRWVSP , BKAX , SPLNAR , SPLNT2 . SP I NTS , PREIIX. 
.  CU.CQ2 

HSTAT  INCLUDE  HSTAT , SPINTG, 5PLVAL. SPFIT , SPLNT2 , CUSC02 , HOWRAL ,  VUNITJ.CtJNI 
, VT , PDER.PAOO , RSOLVE  ,  SPINT2,  SPPLV2  .HORMTS .  ROCOKP .  PfAPY .  PVAl.PINT ,TRIR,  SPL 
.  NFT 

HSTOUT  INCLUDE  HSTOUT 
RESVAV  INCLUDE  REGVAV 
HYDCAL  INCLUDE  HVOCAL 

HY020  INCLUDE  MVD20.  TNOOPT .  GRWLOO,  GRNFRC,  ALAQ.  EXPIWT  .CPFIT  ,«TPCIN, COCOMP 
. , CSOLVE . ATAN3 

T 30 AMD  INCLUDE  T3CAMD , RENT 20 , T IDAHO , CPF IT , ANUPRN , SPF I T , SPINTC , CPLVAL , SPL 
,VAL 

COFOUT  INCLUDE  COFOUT . FtNTSP ,AH3. ROPELN. EXFOR . CPFIT.CPINTO. CPLVAL 
ROSASE  INCLUDE  RDBASE . ROPRIM.NAVNAK.HLLIFT .ROL IFT. SKLIFT.BKLIFT.FNLIFT, S 
. KNF  nc . RDEDOV .  HLEOOV . 8KEDDY . FNEOOY . SKFRSP . EONKSP . RE VAL , CE VAL , SECT , TANAKA 
.  ,  YISC.  SERAB.SERD,  SERE,  FTWO,  F1G5S .  F  IG7  .FIGB.  FIG10.  FIQt  1  .CALRGP.BIL&EK'M 
. AX , Cft INR . F I NTSP . SB  EDDY , SBL I  FT , SPF I T . CPLVAL . CPF  IT 

EOWOTN  INCLUDE  EONOTN, LINIT. SOLVE .CLIP. TRNLAT. RUEVAL .RVSLAT.LSC3F .REVAL . 

, INERST. CSOLVE. CDCOMP.EDMKSP. SKFRSP  FINTSP.ACTF IN 

IRG5EA  INCLUDE  !RGSCA . RMSTOE .KEOEFN.RAOPHS, PRAO. ACRES . ATANSD.ORAO . VELACC 
. . RAOPNA . RELNOT . RHS . TOE . PSPSC . ALGRNQ . SC82 . XMSSC . PSPLC . INTRPL . TEPE AK . 

. FNflAO . LOAD . CPFI T . CPLVAL . CPINTG 
OUTPUT  INCLUDE  OUTPUT . RSTITL 

RAOOUT  INCLUDE  RAOOUT. HLITR.ORGRAO.TFMf  IT  .RAOPNA.  ATAN2D.CPFIT.CPLVAL.SPF 
.IT, SPLVAL. ALGRNQ 

LRAOUT  INCLUDE  LRAOUT.  LRAO.  RAOPHA .ATAN20. CPFIT  .CPLVAL. CPINTO 

RNSOUT  INCLUDE  RMSOUT .  RL l TER ,  FETCH .  $P*.VAL ,  SPFIT , DKMSLW ,  RLITR,  SETSEV , 

.SEVMUT 

HVOCAL  GLOBAL  STELKM 

GLOBAL  OATINP.  lO.PHVSLO.  ENVI0R.RE5PN, STATE, OkON,  APPtlOO.CHW).  INDEX 
GLCTAL  FINCON, LOADS 
EHO 
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APPENDIX  A 

LISTING  OF  UPDATES  TO  SHP  SOURCE  CODE 


•ID  SMP*4 
•0  DAT  IMP. 3 

2  '.RAOPR , QMNOM . KO ,  STATN(  28 ) . NSQFST ( 38 > . 

«D  10.3 

3  ISCARO.BLKFIL.SCRFIL.SPLFIL.LCOFIL.LRAFIL.SEVFIL 
•0  10.5 

a  ISCARO.BLKF IL, SCRFIL, SPLFIL. LCOF IL .LRAFlL.SEVFtL 
•0  INDEX. 2,4 

COMMON  /INDEX/  PFIDX ,  LPFIDX .  RMIOX  .LRKIOX,  SVIDX,  LSVIJX 

INTEGER  LPFIDX. LRNIOX, LSVIDX 

REAL  PF1DX(235)  ,RNIDX(  183).SVIDX(3) 

•0  RESPN.  2 , 3 

COMMON  /RESPN/  NRE5P. IPO!NT( 112 ) . IMOTN( 143  f , ITYPE (  IBS). 

a  ILIN( 183) , ISYM( 1(2) 

•0  R0UE0.3 

2  BKT<  23 ) ,PKS( 28) , S5TR( 2S ) 

•0  SMPB1.3 

a  TAPE3. TAPE4, TAPE  10. TAPE  11. TAPE  12. TAPI13.TAI»t14.TAP*1B) 

•I  PRELIM. 19 

LCOPIL  -  4 
LRAFIL  •  10 
•I  PRELIM. 22 

SEVFIL  ■  14 
•I  READ. 11 
•CA  LOADS 
•CA  FINCON 
•0  READ. 43. 48 

READ  ( I  CARO . 1028 )  OPTN . VLACPR . RAOPR . RLOMPR . LRAOPtt 
WIIITF  i  TPPTN  1.330)  OPTN.  VLACPR  .  RAOPR. RinMM  .LBAOM 

1025  FORMAT  (1615) 

•0  RCA0.4T 

2  4X.  DHVLACPR .  3X .  SNRAOPR  .  4X ,  6HRLOMPR . 4X . 6HL R AOPR/S 1 10) 

•I  RE AO. 63 

C  SPEED  DEFINITION 

IF  (PONITS(I)  ,NE.  METER)  VKMETR  -  VKM61 R/FTMEYB 

METRVK  -  1./VXMETR 

CON  -  VKMETR/S3RT(0RAV*LPP) 

IF  (VKINC  .EO.  0.)  VKINC  -  S. 

IV  -  0 

8  IV  »  IV  *  1 

VK(IV)  •  (  IV-  1  )«Vt1INC 
VFS(IV)  *  VKMETR*VK( IV) 

FRNUM(IV)  -  CON*VK( IV) 

IF  (VK(IV)  , LT .  VKOES  .AND.  IV  .LT.  6)  00  TO  f 
NVK  •  IV 

rHOES  ■  COpj^VKvEo 
«0  READ. 67 

1060  FORMAT  (6F10.4) 

*D  READ  .73,74 

READ  (ICFR0.1020)  NSTATN.NLOADS 
TRITE  (IP) IN. 1100)  NSTATN.NLOADS 
*0  READ. 77. 78 

2  7M0FFSETS//19H  M3.  OF  STATIONS  • , 1 3 ,4X. 6HMLOAOS  ",I3// 

2  3X.7HSTATION.3X. "HMLEVF , 6X , 4HBEAM, SX , 8H0RAFT . 4X , 6HSECARE . 

2  5K.KC0LVL/) 

•l  READ. 136 

IF  (NLOAQS  .EO.  0)  00  TO  BS 
C  READ  WEIGHT  CURVE 

READ  ( ICARO, 1060)  (SVQHT(K).X-I.NSTATN) 

WRITE  (IPRIN. 1333) 

1333  FORMAT  ( //7X ,  12HVE IQHT  CURVE//3X .  7KSTATI0N,  4X .  6HWJI0WT/) 


00  197  K»1 .N3TATN 

197  cowrrwue(,RIM'  10*°*  statn<k),swbhtoo 

c  RE*0  LOCATIONS  (STATIONS)  WHERE  LOAOS  ARE  TO  KF  ril 

Sra  KSffiSS  (XL05TN(K,-K-^0«S)  cut4Tw 

1339  FORMAT  (//4X. 14HL0A0  STATIONS-/) 

as  continue^****  10*°^  <xu>stn(k>.k-i.nu>aos> 

•0  REAO.343.345 

2*JS.  <JCAW5>l°SO)  NFNSET ,  IACTFN,  IFCLCS 
1330  ,a30)  HfNSer.  IACTFN,  IFCLCS 

U30_  *?H*J*T  (  / / /3SH  OATA  CARO  SET  11  -  FIN  PARTICULARS// 
RE20^24|HN,:MSeT,4,ll8HI*CTFW‘4X-CHIFCI-CS/3I10) 

IF  (IACTFN  .60.  0)  GO  TO  132 

!!££*£'  ,0*°)  (FOAIN(IV).IV.I.NW) 

__  „  WITt  (IRRIN.3010)  (VK(IV),IV>1.NVK) 

2010  FO"«AT  (  /23H  SHIP  SPEED  (KNOTS)  o.RFlO.3) 

3030  rn!lfiJlP*,":.30™!  (^INUVJ.IV-I.NWK) 


-.3F10.3) 


».3F  10.3) 


2020  FORMAT  (  32H  FIN  GAIN  FACTORS  -,BF10.3) 

READ  (ICARO.IOAO)  fk  ' 

_ _ WRITE  (IPR1N.3030)  FK 

203°  rowar  J//W,  controller  coeff.  ..jfio  3) 

READ  (ICAR0.1080)  FA  ' 

write  (IPRIN, 30401  F4 

1A4PI  ronu. .  <  » _ _ 

-  »  '**»  aEKWO  COEFFICIENTS  -.SPIO-S) 

READ  (ICAR0.10R0)  FB  ' 

_  "RITE  (IPRIN.3030)  FB 

a°*°  <  /aaH  COMPENSATION  COEFF. •  3F SO. 1 > 

132  If  (IFCLCS  .60.  0)  00  TO  13G 
WRITE  ( IPRIN, 3060) 

3060  FORMAT  (//34X.3CHCORRECTE0  FIN  LIFT  CURVE  SLOPE  I 

S.‘S  IKSiffJ?} 

00  134  K*  1  .NFNSf 

m  s™  !'55S:i™! 

<3°  cSSJ«k,m  """•«.•»•  -  .. 1.10.11 

138  CONTINUE 
•0  READ . 380 

FNIMAQ(K)  •  1. 

«0  REA0^430^443^^^  ^  °'f  fNlMQ^  ’  »• 

*1  READ .511 

c  1*2“  3  FINS 

*0  READ. 3  18 
C  ROLL 

L  *  1 

*0  READ. 534,935 

C  ROU  iE^i^L*CP,,  QT  O  Oa  ST‘T»«^>  «O.EOVLIN»  00  TO  803 
L  «  L  ♦  1 
*0  RE AO. 931 

_  50*  A*  (,*CTFM  «.  O)  GO  TO  908 
C  FIN  8  FIN  VELOCITV 
Ml  ■  1 

IF  (VLACPR  .or.  O)  Ml  •  2 
OO  90S  I7b1.ni 
L  ■»  L  +  1 
IPOINT(L)  -  o 


-  3  FINS 


A-2 


1 

1 

| 

;e 

N 

1 

*2  .** 

»  "X 

H 

V.' 

• 

1 

1  P 

r* 

I 

R 

*  1 

IMOTN(U)  -  0 

■ 

/ 

ITVPE(L)  •  IT 

v  X- 

ILIN(L)  •  .FALSE. 

1 

\  >N 

ISYM(L)  ■  .TRUE. 

i  m 

SOS  CONTINUE 
eoe  NR ESP  -  L 

5  ^ 

GO  TO  aso 

<j  a 

«l  RE AO. 944 

s  *• 

IF  ( IACTFN  .EO.  0)  GO  TO  929 

% 

C  FIN.  FIN  VELOCITY,  AND  FIN  ACCELERATION 

•N 

00  922  IT-1.3 

** 

L  •  L  ♦  1 

Is  >’ 

r 

IPOINT(L)  -  0 

INOTN(L)  >  9 

ITYPE(L)  •  IT 

*'  ■% 

ILIN(L)  «  .FALSE. 

^  *  * 

ISVN(u)  «  -TRUE. 

/ 

922  CONTINUE 

» 

929  CONTINUE 

•I  READ. 904 

< 

s e 

IF  (NLOAOS  . EO.  0)  00  TO  700 

C  LOADS  AT  SPECIFIED  STATIONS 

DO  020  K«1. NLOAOS 

*  • 

C  1-10  (H. SHEAR) 

V  ^ 

C  1-11  (V. SHEAR) 

*  *  *,’ 

*■*,  .  - 

C  1-12  (T.MON.) 

* 

C  1-13  (V.MON.) 

k\ 

C  1-14  (H. MON. ) 

-• 

DO  0to  1-10. 14 

1 1 

IF  (.NOT.  < I ,EQ. 1 1.0R. I.EQ. 13))  GO  TO  910 

L  •  L  ♦  1 

IPOINT(L)  «  K 

k  . 

INIITNJ  L )  •  I 

k",  % 

ITVPE(L)  -  1 

V  *%. 

ILIN(L)  *  .TRUE. 

‘w  -N 

ISVN(L)  •  .TRUE. 

|  „ 

BIO  CONTINUE 

020  CONTINUE 

9  ft 

700  CONTINUE 

*1  READ. 900 

*- 

IF  (NLOAOS  .GT.  0)  LOADS  ■  .TRUE. 

•I  RE AO .SOI 

,N 

IF  (ISK1P  -EO-  1)  ADORES  «  .FALSE. 

«0  READ. 002 

*  '  » 

LRNIOX  «  103 

k  * 

tj  RE AO. 003 

R  f*. 

LSVtOX  *  3 

CALL  OPENMS  ( SEVF IL . SVIOX. LSVIDX .0) 

•t  .*• 

•I  HSTAT . 22 

* 

•CA  LOADS 

•I  HSTAT.20 

; 

REAL  METER 

DATA  METER  /4HNETE/ 

»  *- 

•I  HSTAT.10E 

1 

IF  (NLOAOS  .EO.  0)  GO  TO  09 

V 

C  OBTAIN  LOCATIONS  FOR  LOAD  CALCULATIONS 

V 

DO  09  IP- 1. NLOADS 

<  \*» 

XLS  •  XLDSTN(IF) 

■>*  • 

N1  ■  NSTATN  -  1 

DO  03  K-1.N1 

*■ 

IF  (.NOT.  (XLS.QE.STATN(K)  .t  ».  XLS . LT . STATM(K+1 ) ) )  GO  TO  03 

£  r 

XLOJTN(IP)  •  O.S-(STATN(K)  ♦  STATN(K* 1 ) ) 

X  v 

A~3 

>  V 

r» 

1  sc 

ft  ^ 

: 

1 

GO  70  64 

63  CONTINUE 

64  XLOXPT ( IP  )  •  'C8  -  XLOSTNI IP J»LPP/20 
LSTATN(IP)  ■  NSTATN  +  1  -  K 

69  CONTINUE 

C  COMPUTE  SECTION  MASS 
L  ■  NSTATN  »•  1 
00  66  K“ 1 .NSTATN 
L  -  L  -  1 

IF  (PUNITS(l)  .£0.  METER)  SMASS(L)  ■  SWGHT(K)»  (000 
IF  (PUNITS(I)  .NE.  METER)  SMASS(L)  *  SWQHT(K ) • 2240/08 Ay 
66  CONTINUE 
69  CONTINUE 
•I  HVDCAL.3 
•CA  STELEM 
30  T20AMD.2 

SUBROUTINE  T 20 AMD  (K.PHI20.T20.T30) 

•0  T 20 AMO ■ 12 

COMPLEX  PHI 20 ( 10. 10.4 ) .CTEMP . T20( 10. 1O).T30( 10. 10) 

*0  T20AM0.30 

CTEMP  ■  (O. .O. ) 

•D  T20AM0.46 

CTEMP  ■  CTEMP  *  4TDL(K,K)*N0RM(IN.M.K)*PHI20(ISIGMA,M.UP) 
•O  T20AMO.4R 

TS0{  ISISMA.L)  -  2.0*H-KHu*SI<i»iAiI»»oNAi»xrcTll*CTtMP 
•D  T20AM0 . 5 1 

T30( I  SIGMA , L )  •  T30( I SIGMA , L )  ♦  WTLI «T2D( I3IQMA , L ) 

•I  T 3D AMO . 9 
•CA  STELEM 
•I  T30AMO.14 

COMPLEX  T20( 10.10) 

*0  T3DAM0.39 

CALL  T 20 AMO  (K . PHI20 . T20, T30 ) 

M  •  (K- 1 )• 10 
00  25  L-LMIN.IMAX 
M  •  M  ♦  1 

CALL  CPFIV  (SIGMA. T20{ 1 . L ) . STELEMI 1 . 1 .M) , NSIGMA ) 

29  COWTINUF 
•I  COFCUT.T 
•CA  GEOM 
•I  C0FOUT.S 
•CA  STATE 
•I  COFOUT. 12 
•CA  STELEM 
•I  COFOUT. 14 

COMPLEX  STV(3.3),C0UM(3.3),SF3(25).5H3(25) 

0 1  MENS  ION  SA33( 25) . SB33( 25 ) 

•l  COFOUT. 31 

IF  (OMEGAE  . LT .  SIGKA(O)  OMEGAE  •  SIGMA(I) 

WE  •  OMEGAE 
WE2  -  WE*WE 
♦I  COFOUT. 32 

00  SO  K>1. NSTATN 
SA33(K)  •  O. 

SB33(K)  •  O. 

NPT  -  MOFSET(K) 

IF  (NPT  .LT.  2)  00  TO  60 
M  -  (K- 1 ) *  10  ♦  1 

CALL  AMO  (OMEGAE, STELEM(  1 .  1  ,M)  .  STV.CCHJM) 

SA33(K)  -  REAL(STV(2,2))/( *WC2 ) 

S8330O  •  AIMAO(STV(2.2))/WE 
SO  CONTINUE 


*0  CQF0UT.33 

CALL  AMD  (OMEOAE, TELEM. TV, TL) 

•0  CQFOUT .  30 

CALL  EXFOR  (OMCAA(  IH) .OMEOAE ,  EXCV,  EXCL,  IkJV.HUL ,  H7 , 3F3.SM3) 

*1  COFOUT.37 

IF  ( LOACS  )  ‘-RITE  (LCOFIL)  (SF3(  I  ) .  SH3(  I  ) .  SA33< I  ) . S63J(i  ) ,  !•  t , 

2  NSTATN) 

*0  AW.  3 

flttROUTlNE  AMO  (OMEOAE . YELEM, TV. TL ) 

>0  AMD. II 

COMPLEX  TELEM(4,0, 10) 

•0  EXF0R.2 

SUBROUTINE  EXFOR(QK€OA , OMEOAE . FXV. FXL , HJV, HJL . H7 . F3 ,H3) 

•I  EXFOR. 16 

COMPLEX  CEP, F3( 26 ) ,H3( 29 ) , TF3 , TH3 
•I  CXF0R.22 

ARQLI  ■  '  *N*COSMU 

IF  (ASS(AROLI)  .LE.  TEST)  ARQLI  ■  0. 

•0  EXF0R.33 
•I  EXF0A.34 

IF  (.NOT.  LOADS)  00  TO  210 
F3(K)-(0. .0. ) 

H3(K)-(0. .0. ) 

210  CONTINUE 
•I  EXFOR. 35 

IF  (NOT.  LOADS)  SO  TO  230 
TF3“£KZ*M0RM3«CARQ 
TH3-EK2* 1 1 *TOO«PHI 20( 3 ) 

F3{K)-F3(K)*MTDL( J,K)»TF3 
H3  ( X ) *H3 ( K ) +MTOL ( J , K ) • TH3 
220  CONTINUE 
«D  EXFOR. 67, 69 
*0  EXFQR . 06,90 
•I  EXFQR. 94 

IF  (.NOT.  LOADS)  GO  TO  230 

C  SECTIONAL  FROUOE-KRILOFF  •FORCE",  F3  W/O  CEXP( *1 I“X*X»COS(MU) ) 
F3(K)-2-*CRAVF3(K) 

C  SECTIONAL  DIFFRACTION  "FORCE*  .  H3  M/0  CEXP( -II»K»X*COS(MU)) 
H3(K|«2*W»H3(K) 

230  CONTINUE 
•D  EXFOR. 1O0. 101 
•O  EXFOP. 109 
*0  EXFOR. 114. IIS 
•O  EXFOR. 133 
•  I  C0M0TN.12 
*CA  FINCON 

COMMON  /HULL/  A2S 
•I  EOMOTN.1B 

COMPLEX  ZERO, T*F( 3 ) .CTEMP 
niwtKirai  T«4T(S) 

•I  E0M0TN.20' 

ZERO  •  (0. ,0. ) 

«!  EQMQTN. 41 

WE  •  OMEOAE(IW) 

WE2  •  «•« 

A22  •  RCAL(TL(1.1))/(-«2) 

■22  •  A!MAQ(Tl( 1. 1))/ME 
A26V  »  REAL(TL(1.3))/(-NE2) 

A 26  ■  A 26V  •  ( V/WE2)*B22 
*0  EQMOTN .63,72 

CALL  ROCVAL  (  IV.OMCQA(  IV)  ,0ME0AE(  1M>  ,MRANQ.TLQ.EXC'.«,TLaC.EXCLBC 
2  T44T) 


A- 3 


IF  (IACTFN  .EQ.  0)  00  TO  34 
C  ADO  ACTIVE  FIN  COEFFICIENTS 
OMQE  -  ONEQAC(IM) 

0MQE3  -  OMOE*OMOE 

CALL  ACTF1N  ( I V. ZERO , V . ONCE .0NQE2 , TAF ) 

00  32  I>1,3 

TLOCd.2)  ■  TLOCd.2)  ♦  FOAIW  IV)»TAF(  I ) 

32  CONTINUE 
34  CTENP  •  TLOC( 2,2) 

C  ADO  VISCOUS/' BILQEKEEL  EOOT  OANPINO 
(X)  40  IA*1,NRANQ 
TLOC(2.2)  •  CTENP  ♦  1I*T44T(IA) 

CALL  SOLVE  ( 3.TL0C. EXCLQC.NOTLQ.UL. IP. IPRIN) 

*0  F ML I FT. 23 

O  •  FNINACOO 
•I  FNLIFT.23 

CR  -  XRTF  -  XRTA 

CT  *  XTPF  -  XTPA 

XROC  *  XRTF  -  0. 2S*CR 

XTOC  -  XTPF  -  0.2B*CT 

OX  •  XRQC  -  XTOC 

H  «  SORT(OX*OX  ■>  SPAN* SPAN) 

COSLAM  •  SPAN/H 
SECLAM  ■  1. /COSLAM 

C  LAM  •  ACOS ( SPAN/H )  -  QUARTER  CHORD  SHEEP  AMBLE  IN  RAO! AMS 
•O  FNLIFT.40 

LCS  «  1  .R*PI *EaR/(COSLAM*SORT  I ( EAR*SICLAM)t>*2  ♦  4)  ♦  1.«) 
•I  ROLIFT.24 

CR  *  XRTF  -  XRTA 

CT  »  XTPF  -  XTPA 

XROC  -  XRTF  -  0.3S*CR 

XTOC  -  XTPF  -  0-2S*CT 

OX  •  XROC  -  XTOC 

H  -  SQRT(OX*DX  ♦  SPAN* SPAN) 

COSLAM  •  SPAN/H 
SECLAM  •  1. /COSLAM 

C  LAM  •  AC0S( SPAN/H )  •  QUARTER  CHORD  SWEEP  AMPLE  IM  RADIANS 
•0  ROLIFT.41 

LCS  •  1 . E*PI *EAR/(COSLAM*$QRT( 1 EAA*$ICLAM)**2  ♦  4)  ♦  «.«) 
*1  HLLIFT  2 
•CA  DAT INP 
•0  HLLIFT. 13 

MCHORO  •  LPP 
*0  HLLIFT. IT 
SS  «  O 
SP  •  0 

DO  9  L* 1 .NSTATN 

IF  (L  .EO.  1)  OX  *  { X( 2 )  *  X(  1 )  )/3 

IF  (L  .EO.  NSTATN)  OX  •  (X(NSTATN)  -  X<  NSTATN- 1  )  >/3 
IF  (L.QT.1  .AFC.  L.LT.NS  ATM)  OX  •  (X(L+1)  -  X<L-t))/2 
OX  •  AOS (OX ) 

WPT  •  NOF$ET(L> 

IF  (NPT  .LT.  2)  00  TO  0 
T  •  ABSCd.U) 

A  »  T*OX 
SP  •  SP  ♦  A 
S5  "  SS  ♦  X(L)*A 
9  CONTINUE 
XCP  •  5S/5P 
•D  HLLIFT. 18 

ZCP  •  O. 

‘I  LSC0F.3 


T 

/ 


; 


COMMON  /NULL/  A2R 
•1  LSCOF. R 

COMPLEX  VIV.ZERO.CTENP 
•I  LSCOF . 10 

LOGICAL  HULL 
•I  LSCOF . 12 

ZERO  *  (0. .O. ) 

VItt  -  V/( I ! «OMEQAC ) 

V2W2  •  (V/OMEO*E)»*2 
•I  LSCOF. 33 

MULL  •  -FALSE. 

IF  (MCHORO  EO.  LRP)  HULL  «  .TRUE. 

IF  (.NOT.  HULL)  00  TO  92 
C  HULL 

OF 2  •  ZERO 
0F4  •  ZERO 
0F«  •  ZERO 
SP  •  O 

00  42  L- I. NSTATN 

IF  (L  . EO.  1)  OX  -  (X( 2 )  -  *<1))/S 

IF  (L  -EQ.  NSTATN)  OX  «  (X(NSTATN)  -  XOrtTAYM- I ) )/2 
IF  (L.GT.1  .AM).  L.LT. NSTATN)  OX  •  tX(L*«)  -  X<L-t>)/2 
OX  •  ASS (OX ) 

NPT  -  NOFSCT(L) 

IF  (NPT  -LT.  2)  OO  10  42 
T  «  ABS ( Z( 1 . L  )  ) 

22  •  Z(l.L)/2 
A  •  T  »0X 
SR  •  SR  ♦  A 

F2  >  FZ*OMEOA*(SINaAN»SINMU  -  II'COSOaMS* 

2  CEXR(C!<*(I2  '  » i »» ( i_ ) »C3S!®3 » ) 

CTENR  -  F2«SCN0AN«A 

OF 2  •  0F2  «  CTEMR 

OF«  •  OF*  *  X(L) 'CTEMR  «  VtM*CTENR 
42  CONTINUE 

DF3  *  0F2/SR 
OF«  •  OFB/SR 

CS  -  NKRLA/(LRR*«AAM<DRAFT) 

C*  •  AREAMX/(SEAM«0«AFT) 

CR  •  CB/CX 
OO  TO  «2 
*3  CONTINUE 
•I  LSCOF  SR 
«2  CONTINUE 
*C  LSCOF  . RO 

IF  1  MOT.  >AAL L  I  OMR  •  XCR*XCR*OR22  •  V2S3»Q822 
IF  (FRJU.1  MM  «  (CR»LRR/2)»*2  •  M22  *  V»A2G  ♦  Wi»iS*CS22 
•O  LSCOF  T* 

I*  I  NOT  NA.O  IfCLCCI  S)  •  EICLQIS)  *  Of*  ♦  VI5J*»F2 
IF  (Nxi)  EXCIQCOI  •  EXCLOISI  *  OF! 

*©  RUFfAl  2 

SUMOuTIMR  ROEVAL  1 1 V  .ONSOA  .  QMRRAB  ,  MAASR.TLC,  EF2LR.  I  lOC  .EXCLRC, 

3  T 44T  ) 

•O  *04**1  I 
•C*  FINCON 

COMPAO  UOCi.3).  'XCLG<  JJ.TLOCO.  Jl.lXCLOCCI) 

OINEFFSIlQH  T44TIMIAWR) 

•O  4O4V4L.R0 

TEMR  •  FLCS(K) 

IF  ( IFCLCS  .10.  O  TERM  •  FCLCS(IV.K) 

CALL  LSCOF  (0nt«A.0Mf«Ai.2.FSRAN(K).FRNC»O(K).FAMA(K).TCI». 

•O  ROEVAL.R4 


A- 7 


/ 


2  cau  i to  of  <QNiaA,owcAi,2Jf»«{ic)>o«me<r),rA«*A<c>,TO»t 
isrtfi.12 

ocaoM  •  os(KJ)  •  ra«a) 
slope  «  0- 

if  (otr'  si.  fl-)  step*  •  (Ytsoai)  •  n»  on.  )>/«** 

•o  ursn.zo 

C€#OK  *  U9CCU)  ■  WOCU 

flope  ■»  a. 

IF  (KNUR  -GT.  0.)  SLOr«  »  (TR(«U)  ■  tt'CKDJ/DMCM 
•i  ncvAi.se 

CO  100  IA»1,WMJ«S 

*0  20EVAL.  103 
*D  ROEVAL. 122 
•  I  BDEVAL-U4 

T 44BE  -  O. 

<3  ROE VAL . 127 
•I  ROtVAL-128 

T44T ( (A )  ■  T44SFV  ♦  T44EMV  ♦  T44BC 

lOO  CONTINUE 
•O  RNST0E.25 

CIMCNSION  xio(eii) 

•O  RMS  TOE  . 29 

NIO  •  911 
•I  RMST0E.49 

REWiPO  LCOEIL 
•I  RAOPHS. lO 

COMPLEX  573(25. 30). SH3(J5, 3C> 

DIMENSION  SA3312S, 30). 5833(23. 30) 

•I  RAOPHS. 94 

IF  (IP.EO.O  .AND.  IM.E0.1i)  CALL  FNOAO  (  I  V.M1  .MS.N&TIU  1 .  1 .  IA)  , 
2  RA01(1.IA).PMSI(  I.IAhNMOT.NOMEOA.OMEOAC.tPNS) 

«I  RAOPHS. 63 

IF  (.NOT.  (IP.GT.O  .AND.  < 1M. GE. 10. AND . IM. IE . 14  )  > )  00  TO  IOO 
DO  40  IW*1. NOMEG# 

READ  ( LCOFIL )  ( SF3( * . I V ) . iH3( I . Itf) . SA33l I . l») .  S»33( l . IK) . 

5  I-1.NSTATN) 

40  CONTINUE 

CAU  LRAO  (IM.K1.M2.MaTV.SF3.SH3.SA33.SB33,VFS(IV).CCSMU. 

2  OMEGA. OMEGAE  .  IP.OAOI  .  PMjl  ,FMOT  . NOMEGA  ,  IPHS ) 

*;  rstitl.4 
•Ca  loads 

OIMENSICN  L0A0(2.5),LTVPE(3.2). LUN- T( £ , 3 ) 

REAL  L0*0 . LT VPE . LfNI T 
•I  RrTr.L.29 

DAT*  LOAD  /  4H  H.S.4HHEAR.4H  V . 5 , 4HHE AR , 4H  T  .  4H0RS  .  .  4H  V.B, 

2  ■'HENC..4H  H.B.4HFN0./ 

0  TA  LTYPE  /4HE0RC.4HE  ,4H  , 4HMCMC . 4HHT  ,«H  / 

C4TA  LONIT  /»H  <T.4M0NS).4H  . 4H  (M-.4HTCNI.4H)  .  4H  (FT. 

2  4H-T0N.4HS)  / 

*0  RSTITL. 101 

90  ir  l)K  .NC.  t)  GO  TO  1C 
C  ANTI -ROLL  TINS 
CTITL(  1  )  -  **H 
RTl TL( 2 )  •  4H  FIN 
IF  (IT  .EQ.  1)  JT  «  4 
IF  (IT  .CT.  •)  JT  •  IT 
HO  60  1*1,3 

GO  RTVPE  (  J  >  TVPE(I.JT) 

JT  >  IT  *  3 
00  70  1*1.3 

70  RUNIT(I)  ■  UN!T(I.JT) 


IF 

(  I  r 

•  EO 

1  ) 

FMOT 

-  3MMN 

IF 

(IT 

.£0. 

2) 

FM3T 

•  6HEINVEL 

CF 

(IT 

•  EO. 

3) 

FMOE 

•  6HFINACC 

ENCODE  <  1 10 ,3000, PAHS)  Fim.FNOT 
80  IE  (.NOT.  tlP.aT.O.ANO  < IM. SE . 10. AND. IH. LE . 14 ) J)  CO  Tu  tOQ 
C  LOROS 


«JM  ■  IM  '  9 
RTSTL<«  )  cn  LOAO( 

RVITLC?)  ®  LOACH  3,  JM) 

LT  «.  1 

IE  <SM  .DT.  1 1 J  LT  «  2 
MT  »■  lt 

IF  CLT.EO. 3.AN0.(f>UNITS(1).ME  .METE*))  Ml  *  3 

CiO  81  S«1.3 

RrYPEC )  »  LTVPea.LT) 

RUNtTt?)  -  LUNIT(I.MT) 

82  CONTINUE 

IE  (an  .  %c.  ENCODE  OtO.i03t.PAaS)  PT(IP).Ak.»$tN(IP) 

IF  (JM  .£0.  3)  ENCODE  ( 1.0, 3333 . PARS i  PT(  IP  ) .  XI  OS'?  N(IP  i 

If  (JM  .t3.  3)  ENCODE  (  1  lO,  30yi.  PARS)  PT  < IP) .  XLOSTMl  IP  ) 

IF  'JR  EO.  4)  ENCODE  O 10, 3034. PARS )  PT<IP ) . XLOSTN< IP) 

IF  ( JM  .EO.  3)  ENCODE  ( 1 10. 3033. PARS)  F V ( IP ) . XLD5TM( !° ) 

FORMr.T ( 6HHSHEAR , A3 ,  1 1* .  2VKWIRI 2 .  SFL'At*  FORCE  AT  STATION, 
F0RMAT(6HVSHrAtt.A3. I1K.J9HVERT.  SHEAR  FOJCCE  AT  ST AT  SOW 
FORMAT (4HTM0M.  A3.  tS*  .  JSHTORSIONi'-L  MOMENT  A1  STATION 
AS.  13X.29HVERT.  BEND  MO(S. 

A3,  13X.  29F04SRIZ.  fcXNO.  MOM. 


FORMAT (4HVK0M. 
FORMAT (  4HFM0M, 
CONTINUE 


AT  STATION. 
AT  STATION. 


3031 
3013 

3033 

3034 
3033 

«00 

•I  OUTPUT. 3 
•CA  LOADS 
•I  OUTPUT. 4 

IF  (NLOAOS.QT.O  .and.  LRAOPR.QT.O)  CALL  LRAOur 
*r  rrsout.9 

•CA  LOAOS 
•CA  SEVERE 
•0  RMSOUT . lO 

dimension  xio<aii),vioMe2.*!.t»f«aoi 

DIMENSION  l MOOLI 4 ) , LSVRSP( 13 ) , RSPMHCI 3 . 13) 

•D  RMSOUT . IP 

EQUIVALENCE  (IPOINT.VIO).UIO.MRESP) 

■I  RMSOUT. 20 

DATA  LSVRSP  /3. 3 . 3 . 4. t .t.H. 9.8. 9 . C.t .•/ 

DATA  RSPMK  /4HHEAV.  1NE.4HPITC.  1MH.4H  SUA.1KV.4M  HOL.1HL, 
3  AH  VA.1HH.4MP1VA.1MC.4HP1LA.1HC.4HP3VA.1HC.4HPaLA.1HC. 
a  4HPJVA.  1HC.4KP3LA.  1NC.1HP4VA.  1HC.4HP4LA.1HC/ 

•I  RMSOUT. 33 

NSVRSP  -  13 
•D  RMSOUT. 38 

NIO  •  911 


FS.J.9SK) 

F«.a.»wo 

H.  3.03X1 

Ffi.a.w*) 

n.a.wx) 


•I  RMSOUT. 37 

L  •  LENOTH(RMSFIL) 

M  «  (L-t)/3 

C  sa  -  <39  MEANS  RMS*  II  «*S  GFNfRATIO  BV  SMP8 1 

C  M  •  182  MEANS  RMSFIL  HAS  GENERATED  BV  SMPC4 

K  •  1 

DO  7 30  J»1,3 
DO  730  !• 1 ,M 
K  •  K  ♦  1 
TID(I.J)  •  XIO(K) 

700  CONTINUE 

warsp  >  MRSSP 

00  770  15*1 .NSI8HH 
C  I NO  MOST  PROBABLE  PERIOD 

SW  •  SIQHH(IS) 

IF  (PUHITS(I)  .ME.  METER)  SMH  •  SHHcFTRSTD 
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C  STtfHFICAMT  WAVE  MCIOMT  R  ANDES  BELOW  ABC  IN  HETIRS 

c  sea  state  ( 

IP  (SMH  .LS.  0.88)  PER  -  S.P 
C  SEA  STATE  1 

IF  (SUH  OT.O.t*  .AMD.  SNH.LC.1.3R)  PER  ■  0.0 
C  St 4  STATE  3 

IF  (SVM.OT. 1.34  AMD.  SWT.LE.t.73)  PER  •  7.0 
C  SEA  STATE  4 

SF  (SMH.QT.1TJ  .AW).  SWM.LE.3.34)  PER  •  7.0 
C  SEA  STATS  0 

IF  (SVH.Q.  ,3.2<:  .A*fi.  SRH.L1.1.9T)  PER  -  9.0 
C  SEA  state  g 

IF  (SVN.QT.3.97  .AMD.  SUM. If. 9.94)  PER  -  11.0 
C  SEA  STATE  7 

if  <SVH.flT.1l. 34  .AMD.  SUM. Lt. 12.39)  PER  •  1S.0 
c  sea  state  a 

IF  (SHM. ST. 11.39  .AMD.  SUM.LC. «9.T7)  PER  •  19.0 
C  ARCAT !R  TV  AM  SEA  STATE  C 

SF  <  CRH  .ST.  19.77)  PER  •  19.0 

IF  (PER  .LT.  TunOAL(l))  PER  «  TMODALf 1 ) 

<F  (»ER  .ST.  TROOAL(NTWOO))  PER  -  TK90AL(MTM9D) 

I MOOLI  IS)  «  1 
OO  740  LT-T.MTMOO 

3 F  (&RS(PER-1MX>AL(LT))  .LT.  0.0001)  IMOOUtS)  *  LT 
790  CONTINUE 
7TO  CONTINUE 

iiAi’i*  •  0 

IF  (IMOTW(I)  ME.  1)  ISXPSV  •  1 

c  iskpsv  -  o  all  mu  am  -  output  severe  motion  tables 

C  ISKPSV  »  1  ROLL  MOTION  OWL  T  -  SKIP  SEVERE  WOT  I  on  TABLES 

IF  (ISRA’SV  .(0.  1)  SO  TO  910 
NSVRSP  *  »  ♦  JHO’TLOC 
IF  (NSVRSP  ,9T.  13)  WSVASP  •  19 
CALL  SETSEV  | NSVRSP. LSVP3P, 

NRStMO  •  NSVRSP  *  1 
WSMIMD  •  MSIOWM  ♦  1 
•  30  CONTINUE 
«!  RMSOUT.70 

CALL  STIHDX  (SEVFIL  .RSIMrX.NRSINO) 

00  7  I-1.NRSINJ 
•SIM)K(I)  -  O. 

7  CONTINUE 
•*  RMS0UT.T3 
M  -  O 

IF  (ISKPSV  eg  1)  00  TO  19 
OO  ta  is- I, NSVRSP 
IF  (IR  . WE .  ISVRSP.'LR))  00  TO  1R 
JO  •  Lt 
GG  TO  1« 

1C  CONTINUE 
19  CONTINUE 
-l  RRSOUT.BR 

IF  (JR  .CO.  C)  00  TO  31 

CALL  STUCK  (SEVFIL. SVII0X .NSVIW) ) 

00  4  1*1 .MSMIND 
SVIMDM  I  )  •  O. 

■  CONTINUE 
31  CONTIWiC 
•  I  RWOU7.94 

SVMAAX  •  .  3O3*VM0O«L(  IT0)**3 

IF  (PUNITS(t)  10.  METER)  SMIOWX  •  8WMAA*FTMTR 
•0  RRSOUT. 100. IOI 
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•o  MtsouT.i4a.tsa 

00  120  IVn.NVK 
00  1«0  IH* 1 ,NHC AO 
tcmp  •  ntsTBLUH.iTo.iv) 

IF  ( TEMP  AT.  RtSSMIN)  RMSMIN  -  TCMP 
IF  (TEMP  AT.  RMSMAX)  00  TO  110 
RMSMAX  ■  TEMP 

if  (jr  .so.  o)  oo  to  no 

IF  ( JTO  .HE.  IMOOL( IS ) )  00  TO  110 
IF  (STMMCT  .AW.  IH.GT.13)  00  TO  110 
ttXV  *  IV 
MXH  •  IH 
1(0  CONTINUE 
.*0  CONTINUE 

IF  (OR  .EO.  0)  00  TO  ISO 
IF  (!T0  .HE.  I MOOLI  IS))  00  TO  190 
RSVTOE(I)  ■  MXV 
RSVTOF.(a)  •  MAH 

ic  •  a 

DO  130  IV-1.NVK 
DO  130  IH*  t , NHCAC 
IC  ■  IE  ♦  1 

ItSVTOC  (IE)  •  Mt$TSL(  IH,  ITO,  !  V) 

IE  -  IE  *  I 

RSVTOE(IC)  ■  TOCTBLUH.ITO.IV) 

130  CONTINUE 

C  WRITE  TO  SEVERE  NOTION  FILE 

CALL  WRITMS  (SCVF1L.RSVT0C.IC.1S> 

ISO  CONTINUE 
•D  RMSOUT.  1ST 

IF  (IP.OT.O  .AW.  IM.LE.3)  WRITE  (IPRIH.IOA)  (PTMAUCd .  IP) . 

•  I  RMSOUT  .  2Q2 

!F  UP.ST.e  .AMO.  (IM.OR.  IQ  4*0.  IM  AC  .14)1  WRITE  ( IPRtM.  1«V¥). 

a  xldstmjp) 

10T3  FORMAT  (/sax.THSTATION.FS.  1) 

•  I  RMSOUT  .  lOB 

IF  (IP.OT.O  .AND.  (IM  OC  VQ  ..SMC.  1M.LC.nV)  WRITE  (IPMIM.10B3) 
I  (XT  FORMAT  ( /BTX. 14M( FORCE  /  100  )  ) 

IF  (IP.OT.O  AND.  (IM.OC.ia  .AND.  IM.LE.14l)  WRITC  (IPNIN.10BB) 
IMS  FORMAT  (/S4A.  1«H( MOMENT  /  1Q00O) ) 

■  I  RMSOUT. 333  « 

IF  (IP.OT.O  AND.  (IM.BK.10  .AM).  IMAt.11))  TTMUttUH)  - 
3  TEMRl*S(IM)/100 

IF  (IP.OT.O  . AMD .  (IM.OC.ia  -AMP.  IM.Lt.14))  TCMRHBI  IM)  • 
a  TEMMTSI  IH)/10000 
•I  RMSOUT. a«o 

IF  (IP.OT  O  AND  <  IM.  0(10  AMI.  IR.LC.lt))  T(MU5(  IH)  • 
a  TEMHMSI  IH)/1O0 

IF  (IP.OT.O  AND .  (IM.rC.13  .AM).  IM.LI.14))  TIMRMS(IM)  ■ 
a  f liSwji  iH»/ iuuuu 

•t  RMSOUT. 2** 

IF  (JO  (0.  O)  GO  TO  310 

CALL  STINDX  (  SEVF  I L  .851X35  , MISIND ) 

CALI  WRITMS  (SEVFIL.SWIWX.NSWINO. JR) 

310  CONTINUE 
•I  RMSOUT. 373 

Ir  USKFSV  .(0.  1)  00  TO  4  to 
CALL  ST1WX  I  SEVF  It .  SVIQX  .  LSVIDX  ) 

CALL  WRITMS  ( St VF It  .RS1MM  .«)!)!) ,  IC ) 

4 lO  CONI INUL 

•  I  RMSOUT .  373 


A-ll 
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all  (,**P*V  0>  CALL  SEVNOT  (tSJVESP.  LSVRSP.  R9PNMI  .MONO.  IMQOL ) 

•O  CALrQM.94,13 

C  FIND  HIMIMJM  score  ro«  DCAMISE  CALCULATION  IN  '•ILOCK' 

M2  ■  JS  -  1 
LS  •  K2  -  1 

S«8  •  -  Z(LS.K))  /  (T(M3,K)  -  V(LS.K)) 

J  •  US 

DO  130  1-2. *2 
g  .  j  -  1 
<JS1  -  j  -  1 

SLOPE  -  (Z( J.K)  -  Z(JSI.K))  /  (T(j.K)  - 
IF  (SLOPE  .EO.  0.)  00  TO  140 
IF  (Store  ,«T.  $M)  00  to  i«o 
LS*  JS1 
54®  -  SLOPE 
130  CONTINUE 

C  EXTRAPOLATE  SLOPE  TO  CENTERLINE  TO 

C  ( EXCLUDING  SXE9  OFFSETS) 

140  BKT(X )  -  2(L$,K)  -  SM-T(LS.K) 

IF  (SKT(K)  .LT.  2(1. K))  SKT(K)  - 
•0  S1LCEX.3C 

TLOCAL  »  M$(AKT(K) ) 

-O  SECT. IS 

DIMENSION  AA(3.4),AR( 10) 

•I  SECT. 4* 

AII(I)  -  K 
•0  SECT . 44 

C  .M>T.(ST*TN(M).CO.MKSTN(g.I)))  00  TO  20 
C  ScanrCn  rim  MNIU  RADIUS  OF  THE  OILOE  STANTIN8  FR1M  THE  04  T  ERL  IN! 
AMIN  *  AR(NNN) 

L  -  mooes 
OO  19  NN-2.NNM 
L  •  L  -  1 

m  -  a»(l) 

IF  (*  .or.  SHIN)  00  TO  1 
RMIN  -  R 
IS  CONTINUE 
17  BOX  •  RMIN 
ITSX  «  4 
00  TO  21 
•O  VCOEFN. 13 

OO  HO  1*1.94 
•O  WEOCFW.  14 

no  WCVN(K)  -  0.09  ♦  (I-!)*OWE 
•O  VEOEFN.  «« 

OO  120  1-1.21 

•o  vcoefn. ic 

130  WEWN(K)  -  VEVN(S4)  ♦  I*DVf 
•  »F 


SET  LOCAL  OMAFT 
2(1.  X) 


•COMOECX  LOADS 

CtNWON  /LOADS/  NLOADS .  SMMT (  25 ) ,  SCI 
2  LST ATN( 29 ) 


t( 29) . XLOTTK< lO).rCOXFT(JS), 


•COWECK  FINCON 

CO««ON  /FINCON/  I ACTFN.  IFCLCS.FOAIN(O)  .  FK(  3  ) .  FA(  ])  ,  FB(  3  )  . 
2  FCLCS(S,2) 


-GCMDECK  STCLEM 

COMMON  /STELEN/  STEIEB 
COMPLEX  STELEM(4 , S, 200) 

•COWECK  SEVERE 


I 
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COMMON  /SEVERE/  NRSIND .  Its KNOX ,  WSWI NO,  SWINDX . RSV70E  ,RV, BH 
REAL  RSIN0X( 14 ) . SVINDXt 5 ) . RSVTOE( 40? ) 

INTEGER  RV(13),RH( 13) 

•DECK  ACTFIN 

SUBROUTINE  ACTFIN  ( IV, ZERO .V , OMGE .0MGC3 .TAF ) 

•CA  PHVSCO 
•CA  APPEND 
•CA  RLOOK 
•CA  F INC ON 

COMPLEX  TAF<  3  ).F<3C. CTERM,  ZERO 

00  10  1*1,3 
TAF(I)  •  ZERO 
10  CONTINUE 

FGC  •  ( ( FK(  1 ) -0MGE2*FK(  3 )  )*I  I *OMGE »FK(  2 ))/((( FA(1) -PCJ6C2*FACS}  )•> 
J  II  '0MGE*FA<2)>»<<fB<  1  )  2  *FR  (  3  )  )•!  I  •0MSE*FRI2)  )  ) 

00  30  K'l.NFNSET 

XCP  •  FXCPiK I 

ARM  •  •  FMNCK>(K)/« 

THAT  •  FYHAT'K) 

AP  •  PI*RHO*FSPAM(K>*(FNMCKD(K>/3>**3 
TEMP  •  FLCS(K) 

IF  (IFpLCS  tO.  1)  TEMP  ■  FCLCS(IV.K) 

FZ  -  ( RHO/2 ) »F AREA(K )*TEMP 
SINGAM  •  SIN(FGAMHA(K)*UEGJ)AQ) 

CTERM  •  FGC«( ARM*AP»0MGC2- I 1  •0MGE»(ARM*FZ*3*APJ*V*P?»'#«1f  ) 

Ml  •  1 

IF  (FNIMAQOO  .CO.  3)  Ml  •  2 
C  SIN< IRO-GAMMA)-SIN(GAMMA)  FOR  FIN  ON  STBO  SIDE 
00  20  M* 1 .Ml 

TAF (  1  )  ■  TAF ( 1 )  •  SI NOAM* CTERM 
TAF {2)  •  TAF (2 )  ♦  THAT «CTERM 
r/tr(3)  •  taf ( 3 1  -  ;TMaANaxrp«rTFM 
30  CONTINUE 
30  CONTINUE 

RETURN 

END 

•DECK  FNRAO 

SUBROUTINE  FNRAO  (  IV.NL  .NU.MOTL  .MAQ.PHS.MKIT  . KSMEOA. OMEGAS.  I FN$.) 
•CA  PHVSCO 
•CA  FINCTJN 

COMPLEX  FQC. M8TL( ML  r, NOME OA)  .SETA. ROLI 
DIMENSION  OMEGAS ( NOME GA ) , RAOINOMEOA ) , PHS(MOHEQA ) 

00  10  i-nl.mu 

ROLL  •  M0TL(2. 1  I'RAOOEG 

OMGE  •  OMEGAS ( I ) 

0MGE2  •  0MGE*0MGE 

iQC  •  <  i  Fii(  i  i-u«i»Ei»F*(  ji  )•»  I ’OMu* *rii( 2  )  i/J  {  ir At  i  i*&*5«2"rAi  3? 

2  1 1  •OMGE  *FA(2))*((FS(  1  )  -OMGE  2*  FBI  3  )  )♦!  I  •OMGE*FRI  3)  )  I 
SETA  •  FGAIM(IV)«FGC*R01L 

CAiL  RAOPHA  (RCTA , RAQ( I ) , PN$( I ) . RAOOEG. IPMS) 

10  CONTINUE 

RETURN 

E*4> 

•DECK  LRAO 

SUBROUTINE  LRAO  (  1M ,NL ,NU  ,HJTV,  SF3 ,  SH3  ,  SA33.  SR33 .  V . COS'.V. 

2  OMEGA. OMEGAE. IP. RAO. PHS.NMOT.NOMEOA.IPMSl 
•CA  OATINP 
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OCA  "MYSCC 
•SA  OEDM 
oCA  LOADS 

COMPLEX  ROTV(MKOT.MOMfGA),  IM  .ViWE  .HEAVE  .  HEAVEL  ,  HeAACC  ,  PITCH. 

3  p:mt-. PirACC. VESSEL, V£tfACC.2ERO.n««T.8E$TOfJ,EXCIT. CCP. 

2  HYDRO  ,  LOAO  .  SF I <  2ft .  MBtttCU  f  .  SH4  (39 , NOMiGA  ) 

COMPLEX  S  eW>(5S)tCLEftS(4iaBK=XF.SAB33(SB).CCX«.M¥0,CS.U»i 
DIMENSION  »AO(M^&A),PMS(NOtf£OA).OMEQA(NONEaA),OWEaA£(NWE9A) 
DIMENSION  *A33(SS,NC*CaA).SaJ3(2B.MMtCQA) 

REAL  UCTEP 

DATA  METER  /■vt-KETE/ 

7E«0  “  (C. ,0. ) 

X*>  -  XLOXPT(SP) 

KSTAYN  «  LSTATM(tP)  -  1 
WL>5  ■>  WT.TATN  -  NSTAYN  A  t 
V2  »  V»V 

low  «  tcm 

IF  (PIMP  S<1>  .ME.  METER)  CON  ■  3340 

pi-sks  „<«o.««av 
oo  toe  i^n.,Nu 
%  >  CWE04 m 
w x  »  w»k/ocav 

TEST  -  ,009*TPI/i.lW 
AROU  •  -  WO«COS»*J 

IF  (ASiUAR&LX.*  . L£ .  TEST)  ASOL7  •  Q. 

V£  '■  uH£um£  Hi 

V£?  *  WE*W£ 

l we  -  ii«we 

VIWE  ■  V/IWE 
WE  2  •  V/WEJ 
V2VE3  ■  VS/WFI 
WE  A  Vi!  -  ttJTV(I.I) 

HCAVCL  -  I WE 'HEAVE 
A  ACC  ■»  IWE'Kt'AVEL 
PITCH  •  WITV(I.I) 
riVVEL  -  IWEOFITCM 
PL V ACC  o  rJt*PITVEI. 

VERVll  -  MFAVEL  -  XpoPITVEL 
VEttACC  »  ML' A  ACC  -  XP'OTACC 

C  Ifc'iHT  5  A  TERM 

I’ftLSV  •  JERO 
Ml  ■’  *S»ATH  ♦  1 
00  VO  K-WV.NSTAfM 

StfW(K)  »  SMASS(K)«(.iceXCC  -  ~(K?fPITACC) 
tr  u«  .eo.  t2)  stebhk)  «  -  ix*n)-XP)*5TEMp<»s) 

INERT  "  INERT  ♦  STtMPIKJ 
to  CQftlT  lf*UE 

C  E£SfCiWIM8  TEWS 

DO  30  RESTATE!, MSTAVM 
MPT  *  MOFSETJK) 

SHEAR  >  S»¥(WP.  .  K  > 

STcPPfR)  •  SSEAR ■’(HEAVE  -  L(«)».,’iTCM) 

If  (IW  .EO.  H)  STEMPd!)  »  -  STEKP(K) 

If  (I»  . EO.  13)  *TCW(K)  *  (X(K}-XP)«ATEMP(K) 

70  COW  T I  ML  t 

CALL  C»FIT  (X(K*TATH).STE«*(K*TAT*l),E,,eStf,.'R>l; 

CALL  CPSN’O  (XP,  X(M.;TATN),r{NSTATNJ,NPS,  ElEMS  0.  ,RI*fOR> 

sisroR  »  r>«Ki*aesTO« 
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C  EXCITING  TERM 

DO  30  K’KSTATN.NSTATN 
STEMP(K)  -  SF3(K,  1 )  ❖  SH30<.I) 

IF  (IN  .CO.  13)  STEKP(K)  •  -  ( (X(K)-XP)*$TZMP(K)  * 

3  VIWE«*W3(X,I)) 

30  CONTINUE 

CALL  CPF  IT  (X(KSTATN) , STEMP(KSTATN) ,  El  CMS ,NPS) 

CALL  CPINTO  (XP.X(NSTAYN)  .X(KSTATN)  .NPS.ELEMS.AMi.!.  EXCIT) 

IF  (.NOT.  IN. CO. 11)  GO  TO  33 

CALL  CPF  I T  <X(KSTA'N),SH3(KSTATN.!).ELEMS,NPS) 

CALC.  CPLVAL  (X(KSTATN),.<N>S.ELEMS,XP,eXF,CDUN,IELMl 
CEP  •  CEXP(  II*XP-AR0LO 
EXCIT  ■  EXCIT  *  VIME»CEP*EXF 

39  EXCIT  •  RHO*EXCir 

C  HVOROOVNANIC  TERM 

00  40  K-KSTATN.NSTATN 
A33  •  $A33(K.I) 

■33  -  S)?3(K,I) 

SAB33(K>  •  A33  <■  II*B33 

IF  (IM  . £Q.  11)  STEUP(K)  «’  -  <A33*(H£AACC-X(K)*BITKC)  ♦ 
a  B33*(HrAy£;,-X(K)»PITVEl)  -  VN£2»833»FITACC  ♦  V*A13«PITVEL) 

IF  (IM  .EO.  13)  5TEM»(K)  -  ( K(K ) -XI' ) « ( A33«(HCAACC-X(K)*PIYACC) 
a  V  S33«(HEAVEL-X(K)«PtrvEU)  *  (V*A33»VERVEL  ••  WE2«BOT«VltaACC  •. 
a  -  V2WE:»(A3VPITACC  ♦  B3J*PIWEl)) 

40  CONTINUE 

CALL  CPF  IT  (^(KSTATN).STCMP(KSTATN),ELr.NS.NP3) 

CALL  CPINTG  (XP.X(NSTATK).X(X3TATN),M»S.  ELENS.  O..WMO) 

IF  (  A  flf .  IM.  FO  11)  OB  TCI  48 

CALL  CPFIT  (X(KSTATri) „SAS3S(RSTmTM) ,ELEtSS,K?3) 

CALL  CPLVAL  (X(KSTATM) .WP*. ELEKS.Xr-.HYD.CDWt. IfLM) 

A33  -  BEAL(MVO) 

833  -  AIMAQ(HVO) 

HYDRO  «  HVORQ  -  (VA33«VtRV£L  -  WCa*S3J*VCRAC«  - 

a  vavEa«(A33^PiTfcCC  ♦  B33*>pitvel)) 

49  CONTINUE 

CSIXJ  «  RESTOM  *  EXCIT  ♦  HVSMQ 
LCIAO  »  INERT  -  CSUN 
LOAD  ■  LOAO/CON 

CALL  RAOPMA  (  LCAO  ,CAO(  I  )  .  PHS(  I )  „  BAOTEa.  I  PUS  3 
100  CONTINUE 

RETURN 

END 

'DiCK  LftaOuT 

S'JbROUTIME  LRAOUr 

>*ca  cat  imp 

•  CA  OECW< 

•CA  STATE 

•  CA  PNVSCO 
•CA  CNVIOR 
•CA  10 
•CA  LOWS 

CUWLEX  M0TV(3.3O).MOTU3,3O.3).IWV(3.^).NO',<S  2CJ.NT(*UJ. 
a  SF3'  29.30). 5K3(29.  JO) 

DIM?  NS  I  ON  S*33(IJ,30).S*T3(2S,H3}.t.L'VJiA^OO).V‘*TK.«H}iaS, 

2  VSHvS(X»).VEAL<A0(3O).VB»PM513Oy 
REAL  ME. EM 
CAT*  METER  / 4HMETE/ 

IF  (PUNITSd)  EO.  METER)  UNITE.  *  7t*i- TOWS 


A-1S 


IF  (FUNITS(I)  .F'E.  METER)  UNITS  «  7MFT-T0HS 
SEWINt!  LCOFIL 
RFWIM!  OMOFIL 

READ  (UWiFIL)  TITle  .NVX » NOW , NOME  (LA  .ONEGA  .FIRXNO,  RLANQ,  VRT ,  LAT , 
a  A  IV.iRl  S .  LP1» .  B  E  AH ,  OR  A  F  T .  D I 5PLM .  <3N .  0ELGM . K0 . KROLL . LCS . ORA V . RHO , 

2  VKDES.VK INC. DHLWL 

WRIT E  (LRAFIL)  TITLE. NOMEOA, ONEGA. NVK.NMU.LPP, BEAN, DRAFT. OISPLM, 
J  m . OELGM .I<a , KftfllL .  U.B , OSLttl.  .(WAV, NSTATN , STATN, FM.OAOS , SWOHT , SIM  5 1 
*  XLOSTW.XI.OXPT.X 
IKS  300  IV*I.NWK 
DO  tOO  IKol.MMJ 

READ  mm!  VXVtOTS.HtAONU.ONEOAE 
IF  ( VRT )  RE  AO  (ORGFIL)  MtJTV 
IF  (LAV!  READ  <0S3fIL)  MI3TL 
IF  (AL-3SES)  READ  (ORUFIL)  HJV.HJL.H7 
t-OMl  «  <m.  -  MEAiDMO 
CD5MN  ••  Ci)S(MU(IH,  IV)  ) 

DO'  to  (W«f,NDMf<M 

HEAD  (LCOFIL)  (SF3(I,IW),  SH3(  I .  IW) ,  SA33(  I ,  Itf ) .  1B3J( I , IW) , 

£  I«t.NSTAT») 

tO  CON?  KMUL 

PI>  too  IF"  UNLOADS 
IM  *  «t 

CALL  LRU!  (IN,  1 .NONEQA ,MOTV. SF3.  SH3, SA33 .  SB  33 .  VFS(  IV)  .COSMJ, 

3  OMEOA.OHCttAE.IP.VSFRAO.WSFPHS.a.NOMEQA. 1) 

IM  »  13 

CALL  IRAO  (IN.  1.N0MEQA.MQTV.3F3.SH3.S&33.SJS3a,VfS(lV).C03eU. 
a  ONIOA.OMCOAE.  IP,  V6MRAO.VBKPMS.3.NOMEOA.  1  ) 

WHITE  (IPRIN. 1000)  TITLE , XLDSTN( IP) . VXNOTS.HDNG 
10CU  FORMAT  (  •  1  »/2*X  .  20A4/// I3XJL0AD  RESPONSE  AMPLITUDE  OPERATORS* 
a  *  (RADS)  AW  PHASES*///5OX*STATI0N*F».  1///S5X*SHIP  STEED  **, 

2  F3.C->  KN0TS*/S3X*SHIP  HEADING  **F8.0*  DEGREES*) 

ItitlTE  (IPRIN.  1010)  UNITS 

1010  FORMAT  <//20X*V. SHEAR (V3)*9X *V. MUM. (VS)*/2X*0MEQA  0MC0AE*4X. 

2  2(*AHPL.  PHASE*. 4X)/4X*RPS*4X,*RPS*4X*  T0NS*#X*0E«*4X,A7. 

£  3X*OE0*/) 

00  30  IW«1  .M3MEG4 

WRITE  (IPRIN. 1020)  0ME8A( IW) .0ME8AE( IW) . VSFRAO( IW) . VSFPHS( IW) . 

2  VBM«AO(IW).VBMPHS(IW) 

1020  FORMAT  (2F7.3,3(  1PE12.4.0PF7. 1)) 

20  CONTINUE 

WRITE  (IPRIN. 1030) 

1030  rQtfMAT  (//2.V»NQT5:  HSiBSNS  C5NVEWYI0 Ji:  0  B£S=M€£",  SS  5E5=s 
~3  »  STUD  BEAM.  1*0  OCO*  FOLLOWING  SEAS.*) 

WRITE  (LRAFIL)  XLD5TN(  IP)  ,  VKNOTS.FiDMO.OMEGAE  ,  VSFRAO,  VSFPHS. 

2  VBNKAO.VBMPMS 
lOO  CONTINUE 
IOC  CONTINUE 
300  CONTINUE 


REWIND  ORGFIL 


REWIND  LRAFIL 


RETURN 

END 

•DECK  SETSEV 

SUBROUTING  SETSEV  (NSVASP.LSVIISP) 
•CA  RESPN 

DIMENSION  LSVRSP(MSVRSP) 


00  1&3  LR*?,MSVSSP 
DO  140  1 13*1 .NRESP 
IP  »  IPOINT(SR) 


IN  •  IWOTN(IR) 

IT  ■  I  TYPE  (lit) 

00  TO  (10. 29. 30.40. SO. 60. 70. 60. 90. 100, 110. 130. 130). LR 
10  IF  (.NOT.  ( IP.EO.Q  .AM).  IN.E0.3  .AM).  IT.EQ.1))  00  TO  140 
C  HEAVE 

00  TO  190 

20  IF  (.NOT.  (IP.CO.O  -AM).  IW.EO.B  -AND.  IT.E0.O)  00  TO  140 
C  PITCH 

GO  TO  ISO 

30  IF  (.NOT.  ( 2P.EQ.0  .ANO.  IN.EQ.2  .AND.  XT.E0.1))  00  TO  140 
C  SWAY 

GO  TO  1JQ 

40  IF  (..NOT.  ( IP.CO.O  .AM).  IN.CQ.4  .AND.  IT.EQ.1))  00  TO  149 

C  ROLL 

GO  TO  ISO 

SO  IF  (.NOT.  (IP.EQ.O  .AND.  IN.EQ.6  .AND.  IT.EQ.1))  90  TO  140 

C  YAW 

GO  TO  ISO 

60  IF  (.NOT.  (IP.E0.1  .AND.  IM.E0.3  .AMD.  ST.EO.S))  03  TO  140 
C  VERTICAL  ACCELERATION  AT  POINT  1  (PI) 

GO  TO  ISO 

?0  IF  (.NOT.  (IP. CO-1  .AND.  IN.E3.2  .AM).  IT. CO. 3))  QO  TO  140 
C  LATERAL  ACCELERATION  AT  POINT  1  (PI) 

GO  TO  ISO 

SO  IF  (NOT.  (IP.EO.a  .AND.  IN.E0.3  .AND.  1T.E0.3))  UU  VO  140 
C  VERTICAL  ACCELERATION  AT  POINT  2  (P2) 

GO  TO  150 

90  IF  (.NOT.  (IP.EO.a  -ANO.  IN.E0.2  .ANO.  IT.EQ.3))  GO  TO  149 
C  LATERAL  ACCELERATION  AT  POINT  2  (P2) 

GO  TO  190 

100  IF  (.NOT.  (IP.EQ.3  .AND.  IN.EQ.3  .ANO.  IT.E0.3))  QO  TO  140 
C  VERTICAL  ACCELERATION  AT  POINT  J  (P3) 

GO  TO  ISO 

110  IF  (.NOT.  ( IP.E0.3  -AM).  IM.EQ.2  .AND.  IT.C0.3))  flU  TO  140 
C  LATERAL  ACCELERATION  AT  POINT  3  <P3» 

GO  TO  180 

120  IF  (.NOT.  (IP.EO.A  .AM).  IN.E0.3  .AND.  IT.EQ.3))  BO  TO  140 
C  VERTICAL  ACCELERATION  AT  POINT  4  <P4) 

GO  TO  150 

130  IF  (,M) T.  (IP.EO.a  .AND.  IN. CO. 2  .AND.  IT.EO.3))  00  TO  140 
C  LATERAL  ACCELERATION  AT  POINT  4  (P4) 

GO  TO  150 

140  CONTINUE 
130  LSVRSP(LR)  «  IR 
160  CONTINUE 

RETURN 

END 

•DECK  5EVM0T 

SUBROUTINE  SEVMOT  (NSVRSP .  LSVRSP .  RSPMN . H9NQ.  IMOOL ) 

•CA  OATINP 
•CA  INOFX 
•CA  GEOM 
•CA  PHVSCO 
•CA  10 
•CA  ENVIOR 
•CA  SEVERE 

DIMENSION  RSV( 13. 13),T0E( 13. 13),TENV( 13 ) , TENW( 13) , TEN*( 13) , 
a  TEMT(«3).LSVRSP(13).RSPMK(2.13).l«Ma(34).IM00L<4) 

INTEGER  TENT 
REAL  METER 
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OATA  METER  /4HMETE/ 

NHEAO  «  34 
N1  •  NHEAO  ♦  1 
NDATA  -  2  +  N1*NVK»2 
DO  500  IC-1.2 
DO  400  IS*1 .NSIGWH 
LT  -  IMOOL(IS) 

00  340  IR* 1 .NSVRSP 
DO  200  JR* 1. NSVRSP 

CALL  FETCH  ( IC. JR , IS , RSVTOE . 5VI0X.RSIN0X, SWINDX.NDATA. LSVIOX. 

2  NRSINO.NSWINO.SEVFIL) 

IF  (IR  .OT.  I)  GO  TO  10 
RV( JR)  -  RSVTOE ( 1 )  ♦  .001 
RH( JR )  •  RSVTOE (2)*  .001 

10  IF  (JR  . QT .  1)  GO  TO  20 
IV  •  RV(IR) 

IH  •  RH! IR) 

20  IE  ■  3  ♦  (IH-1)«2  ♦  ( IV- 1 )*NHEA0*2 
RSV(JR.IR)  -  RSVTOE(IE) 

TOE(JR.IR)  *  RSVTOE(IEO) 

200  CONTINUE 
300  CONTINUE 

WRITE  (IPRIN.1000)  TITLE 

lOOO  FORMAT  (*1*/3RX.20A4///4RX»S  EVCRC  NOTION  TAIL* 

3  *E*) 

IF  (IC  .EO.  1)  WRITE  (IPRIN.IOIO) 

IF  \  IC  .EG.  2)  WRITE  OPwiN.  1030) 

1013  FORMAT  (//SOX. 1 IHLONOCREJTCO) 

1020  FORMAT  <//60X. 12HSH0RTCRESTED) 

IF  (PUNITS(I)  .ME.  METER)  WRITE  (JPRIH.1030)  SIQWH(IS) 

IF  (PUNITS(I)  .EO.  METER)  WRITE  (IPRIM.1040)  SI0WHOS) 

1030  FORMAT  ( / 42X*SCA  STATE:  SIGNIFICANT  WAVE  KEIOHT  ■>«.!•  FEET  •) 
1040  FORMAT  ( / 42X*S£A  STATE:  SIGNIFICANT  WAVE  HEIGHT  **FS.2*  METERS*) 
WHITE  (IPRIN.IOSO)  TmooAL(LT) 

IF  (NSVRSP  -EO.  5)  GO  TO  SO 
NP  •  NSVRSP  -  3 
NP  '  NP  /  2 
WRITE  (IPRIN. 1025) 

1025  FORMAT  (//54X*P0INT  LOCATIONS:*) 

DO  50  IP*  1 . 14* 

WRITE  (IPRIN. 1026)  IP , (PTNAMEt I , IP) . I- 1 .«) , XPTLOC( IP) . 

2  VPTLOC(IF).2PTL0C( IP) 

102C  FORMAT  (22X*P*I1*-  •«A4.2X.*XFF  *»F7 , 2 , 2X*VCl_  ••PT.2.2X, 

2  *Z8L  **F7.2) 

50  CONTINUE 
60  CONTINUE 

1050  FORMAT  (54X-M0 C'L  WAVE  PERIOD  **F4.0*  SECONDS*) 

WRITE  (  IPRIN,  1C  :*.)  (STATNM(I). 1*1,3) 

1055  FORMAT  (//40X.3>-.«  VALUE  /  ENCOUNTERED  MODAL  PERIOO  (TOE)*) 

WRITE  (IPRIN.  104.)  ((ASPMKd.IR).  1*1, 2).  IR-1,  NSVRSP) 
lORO  FORMAT  ( // ATX ‘MAXIMUM  RESPONSES  AND  CONDITIONS*/ IX , 130( IH-)// 

2  •  RESPONSE  *. 13(4X.A4,A1)) 

00  310  IR*1. NSVRSP 
IV  •  RV(IR) 

EH  -  RH(IR) 

TENV(JR)  «  VK(IV) 

TEFHUR)  ■  HDNO(IH) 

TEMR(IR)  -  RSV(IR.IR) 

IF  (]R  .QT.  5)  TERROR)  •  TERROR)  •  100 
TEMTOR)  •  TOE ( IR.  IR ) 

310  CONTINUE 

WRITE  (IPRIN,  1070)  (TEMR(  IR)  ,  TENT(  IR)  ,  (R*1  .NSVRSP) 


1070  FORMAT  (/*  (MAX-RSV)/TOE*. 13( 1X.FB.2. 1H/, 12) ) 

Wilt  ( IPRIN,  1040)  (TEMV(  lfl)  ,  IR*  1 , NSVR5P ) 

1000  FORMAT  (»  AT  SFEED  (KNOTS)*. F«. 1 , 12F9. 1 ) 

WRITE  (IPRIN. 10*0)  ( TEMH( IR ) , IR* 1 . NSVRSP ) 

10*0  FORMAT  (•  AT  HEAOINO  (OEO) • , IR.O, 13FP.O) 

WRITE  (IPRIN. 1100)  ( ( RSPNME (I,JR),I*1,2).J9*1, NSVRSP ) 

1100  FORMAT  (//S4X* ASSOCIATED  RESPONSES*/ IX, 130( 1H- )// 

2  •  MAX.  SPEED  /*/■  RESPN.  HI AO IMS* . 3X,A4,A1,12(4X,A4,A1)) 

WRITE  (IPRIN, 11 10) 

1110  FORMAT  (IX) 

00  330  IR* 1 .NSVRSP 
IV  •  RV(IR) 

IM  *  RH(IR) 

MV  ■  VX(IV)  ♦  .001 

IM  -  HDNO(IH)  ♦  .001 

IF  (IR.EO.e  .OR.  IR.EO.fi  .OR.  1R.EQ.10  .OR.  IP. (0.12) 

2  WRITE  (IPRIN. 11 10) 

DO  320  JR* 1. NSVRSP 
TEMR(JR)  *  RSV(JR.IR) 

IF  (JR  .QT.  9)  TEMO(ja)  -  TEMR(Jft)  •  100 
TCMT(JR)  •  TOS(JR.IR) 

320  CONTINUE 

WRITE  (IPRIN. 1120)  (RSPNME (I.IR).I*1,2) ,MV . MH, (TEMR( JR) ,TEMV( JR) . 
2  JR* 1. NSVRSP) 

1120  FORMAT  ( IX, A4. A1 . 2X . 12 . 1H/ . 13. 13(F*. 2 . 1H/ . 12) ) 

330  CONTINUE 

WRITE  (IPRIN. 1130) 

1130  FORMAT  (//2X*N0TES:  I)  RESPONSES  ARE  IN  PHYSICAL  UNITS:*/ 

2  22X.*HE«VE  AND  SWAY  ARE  IN  WAVE  NEIOHT  UNITS;  PITCH.  • 

2  'ROLL,  AM)  YAW  ARE  IN  OCQREES ;  */22X ,  *AN0  THE  POINT  VERTICAL  • 

2  *AND  LATERAL  ACCELERATIONS  ARE  IN  UNITS  OF  Q-S  *.1H*.*  ICO.*) 
WRITE  (IPRIN. 1140! 

1140  FORMAT  («X*2)  POINT  LOCATIONS:  XFP  IS  IN  STATION  MMERS;  • 

2  *VCL  AM)  ZB l  ARE  IN  WAVE  HCIQNT  UNITS.*) 

WRITE  (IPRIN. 11MQ) 

1180  FORMAT  («X*3)  HEAOINO  CONVENTION:  O  OEO*HEAO„  *0  DfO-STWO  REAM. • 

2  »  1BO  OEO«FOLLOMINO  SEAS.*) 

400  CONTINUE 
BOO  CONTINUE 

RETURN 

EM) 
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OTNSROC  ISSUES  THREE  TYPES  OF  REPORTS 

1.  OTNSROC  REPORTS,  A  FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH¬ 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 


2.  DEPARTMENTAL  REPORTS.  A  SEMI  FORMAL  SERIES.  CONTAIN  INFORMATION  OF  A  PRELIM¬ 
INARY,  TEMPORARY.  OR  PROPRIETARY  NATURE  OR  OP  LIMITED  INTEREST  OR  SIGNIFICANCE. 

THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 


X  TECHNICAL  MEMORANDA.  AN  INFORMAL  SERIES.  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  ANO  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN¬ 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AMD  THE 
NUMERICAL  COOE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSNDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE4V-CASE 
BASIS. 


DEPARTMENT  OF  THE  NAVY 

NAVAL  SURFACE  WARFARE  CENTER 
CARDCROCK  DIVISION 


CARMRCCK  DIVISION  HEAOCKJARTERS 
DAVID  TAYLOR  MOOCL  SASIN 
BCTHESDA,  MO  200S4-5000 


Defense  Technical  Information  Center 
Bldg  #5,  Cameron  Station 
Alexandria,  VA  22304-6145 


IN  REPLY  REFER  TO: 

5605 

504 

2  9  SEP  1995 


Gentlemen: 

Carderock  Division,  Naval  Surface  Warfare  Center,  report  DTNSRDC  SPD- 
0936-04,  "SMP84:  Improvements  to  Capability  and  Prediction  Accuracy  of  the 
Standard  Ship  Motion  Program  SMP81 by  William  G.  Meyers  and  A.  Erich 
Baitis  has  been  Approved  for  Public  Release  as  of  1 9  June  l^Q^PJaasa^-o^ 
amend  your  copies  of  the  report  to  reflect  this  changev^i^^50^15®r 

Sincerely, 


<24.  A.  Ufa-' 

W.  B.  MORGAN 

Head,  Hydromechantee  Directorate 
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DEPARTMENT  OF  THE  NAVY 

NAVAL  SURFACE  WARFARE  CENTER 
-  CARDEROCK  DIVISION 

/?/> .  S0?6J3'7? 

Defense  Technical  Information  Center 
Bldg  #5,  Cameron  Station 
Alexandria,  VA  22304-6145 

Gentlemen: 


CARDEROCK  DIVISION  HEADQUARTERS 
DAVID  TAYLOR  MODEL  BASIN 
BETHESDA,  MD  20084-5000 


5605 

504 

2  9  SEP  1995 


Carderock  Division,  Naval  Surface  Warfare  Center,  report  DTNSRDC  SPD- 
0936-04,  "SMP84:  improvements  to  Capability  and  Prediction  Accuracy  of  the 
Standard  Ship  Motion  Program  SMP81,"  by  William  G.  Meyers  and  A.  Erich 
Baitis  has  been  Approved  for  Public  Release  as  of  19  June  1990.  Please 
amend  your  copies  of  the  report  to  reflect  this  change.  .(nt>-ko<XfOi9) 

Sincerely, 


(bt.  A- 


W.  B.  MORGAN 

Head,  Hydromechanics  Directorate 


